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TITLE OF THE INVENTION 
BIOREACTIVE ALLOSTERIC POLYNUCLEOTIDES 

CROSS-REFERENCES TO RELATED APPLICATIONS 

This application is a continuation-in-part of co-pending U.S. application 
serial number 09/331,809, filed June 18, 1999 as a national phase entry of 
PCT/US97/24158, filed internationally December 18, 1997 and claiming priority 
5 benefit of U.S. application serial number 60/033,684, filed December 19, 1996 and 
serial number 60/055,039, filed August 8, 1997. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH 
The invention was made with partial government support under NIH grant 
GM59343. The government has certain rights in the invention. 

10 BACKGROUND OF THE INVENTION 

Field of the Invention . This invention relates primarily to functional DNA 
polynucleotides that exhibit allosteric properties, and to catalytic RNA and DNA 
polynucleotides that have catalytic properties with rates that can be controlled by a 
chemical effector, a physical signal, or combinations thereof. Bioreactive allosteric 
15 polynucleotides of the invention are useful in a variety of applications, particularly as 
biosensors. 
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Biosensors are widely used in medicine, veterinary medicine, industry, and 
environmental science, especially for diagnostic purposes. Biosensors are typically 
composed of a biological compound (sensor material) that is coupled to a transducer, 
in order to produce a quantitative readout of the agent or conditions under analysis. 
5 Usually, the biological component of the biosensor is a macromolecule, often subject 
to a conformational change upon recognition and binding of its corresponding ligand. 
In nature, this effect may immediately initiate a signal process (e.g., ion channel 
function in nerve cells). Included in the group of 'affinity sensors' are lectins, 
antibodies, receptors, and oligonucleotides. In biosensors, ligand binding to the 
10 affinity sensor is detected by optoelectronic devices, potentiometric electrodes, field 
effect transistors (FETs), or the like. 

Alternatively, the specificity and catalytic power of proteins have been 
harnessed to create biosensors that operate via enzyme function. Likewise, proteins 
have been used as immobilized catalysts for various industrial applications. The 

15 catalytic activity of purified enzymes or even whole organelles, microorganisms or 
tissues can be monitored by potentiometric or amperometric electrodes, FETs, or 
thermistors. The majority of biosensors that are commercially available are based on 
enzymes, of which the oxidoreductases and lyases are of great interest. It is nearly 
exclusively the reactants of the reactions catalyzed by these enzymes for which 

20 transducers are available. These transducers include potentiometric electrodes, FETs, 
pH- and 0 2 -sensitive probes, and amperometric electrodes for H 2 0 2 and redox media- 
tors. For example, the oxidoreductases, a group of enzymes that catalyze the transfer 
of redox equivalents, can be monitored by detectors that are sensitive to H 2 0 2 or 0 2 
concentrations. 

25 Enzymes are well-suited for application in sensing devices. The binding 

constants for many enzymes and receptors can be extremely low {e.g., avidin; K d = 10" 
15 M) and the catalytic rates are on the order of a few thousand per second, but can 
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reach 600,000 sec 1 (carboanhydrase) (45). Enzymes can be monitored as biosensors 
via their ability to convert substrate to product, and also be the ability of certain 
analytes to act as inhibitors of catalytic function. 

Organic chemistry and biochemistry have reached a state of proficiency 
5 where new molecules can be made to simulate the function of protein receptors and 
enzymes. Macrocyclic rings, polymers for imprinting, and self-assembling monolayers 
are now being intensively investigated for their potential application in biosensors. In 
addition, the immune system of animals can be harnessed to create new ligand-binding 
proteins that can act as artificial biorecognition systems. Antibodies that have been 
10 made to bind transition-state analogues can also catalyze chemical reactions, thereby 
functioning as novel 'artificial enzymes' (36). The latter examples are an important 
route to the creation of biosensors that can be used to detect non-natural compounds, 
or that function under non-physiologic conditions. 

Description of Related Art . In nature, RNA not only serves as compo- 
15 nents of the information transfer process, but also performs tasks that are typically 

accomplished by proteins, including molecular recognition and catalysis. A seemingly 
endless variety of aptamers, and even DNA aptamers can be created in vitro that bind 
various ligands with great affinity and specificity (17). Nucleic acids likely have an 
extensive and as yet untapped ability adopt specific conformations that can bind 
20 ligands and also to catalyze chemical transformations (16). The engineering of new 
RNA and DNA receptors and catalysts is primarily achieved via in vitro selection, a 
method by which trillions of different oligonucleotide sequences are screened for 
molecules that display the desired function. This method consists of repeated rounds of 
selection and amplification in a manner that simulates Darwinian evolution, but with 
25 molecules and not with living organisms (4). One drawback to the use of existing 
enzymes as biosensors is that one is limited to developing a sensor based on the 
properties of existing enzymes or receptors. A significant advantage can be gained if 
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one could 'tailor-make' the sensor for a particular application. It would be desirable 
to employ nucleic acids to create entirely new biosensors that have properties and 
specificities that span beyond the range of capabilities of current biosensors. 

In vitro selection has been the main vehicle for new ribozyme discoveries in 
5 recent years. The catalytic repertoire of ribozymes includes RNA and DNA phospho- 
ester hydrolysis and transesterification, RNA ligation, RNA phosphorylation, alkyla- 
tion, amide and ester bond formation, and amide cleavage reactions. Recent evidence 
has shown that biocatalysis is not solely the realm of RNA and proteins. DNA has 
been shown to form catalytic structures that efficiently cleave RNA (5,7), that ligate 
10 DNA (10), and that catalyze the metallation of porphyrin rings (24). As described 
herein, self-cleaving DNAs have been isolated from a random- sequence pool of 
molecules that operate via a redox mechanism, making possible the use of an artificial 
DNA enzyme in place of oxidoreductase enzymes in biosensors. In addition, these 
DNA enzymes or 'deoxyribozymes' are considerably more stable that either RNA or 
15 protein enzymes - an attractive feature for the sensor component of a biosensor device. 



BRIEF SUMMARY OF THE INVENTION 

It is an object of the invention to provide examples of RNA and DNA 
sensing elements for use in biosensors, including polynucleotides attached to a solid 
support. Both RNA and DNA can be designed to bind a variety of ligands with 
20 considerable specificity and affinity. In addition, both RNA and DNA can be made to 
catalyze chemical transformations under user-defined conditions. A combination of 
rational design and combinatorial methods has been used to create prototype biosensors 
based on RNA and DNA. 
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These and other objects of the invention are accomplished by the present 
invention, which provides purified functional DNA polynucleotides that exhibit 
allosteric properties that modify a function or configuration of the polynucleotide with 
a chemical effector, a physical signal, or combinations thereof. The invention further 
5 provides purified functional polynucleotides having catalytic properties with rates that 
can be controlled by a chemical effector, a physical signal, or combinations thereof. 
Some embodiments are enzymes exhibiting allosteric properties that modify the rate of 
catalysis of the enzyme. In addition, the invention encompasses biosensors comprising 
bioreactive allosteric polynucleotides described herein. 

Examples of chemical effectors include, but are not limited to, organic 
compounds such as amino acids, amino acid derivatives, peptides, nucleosides, 
nucleosides, nucleotides, steroids, and mixtures of organic compounds and metal ions. 
In some embodiments, the effectors are microbial or cellular metabolites or compo- 
nents of bodily fluids such as blood and urine obtained from biological samples. In 
other embodiments, the effectors are pharmaceuticals, pesticides, herbicides, and food 
toxins. Physical signals include, but are not limited to, radiation and temperature 
changes. 

The invention also provides methods for determining the presence or 
absence of compounds, or compound concentrations in biological, industrial, and 
20 environmental samples, and physical changes in such samples using bioreactive 
allosteric polynucleotides of the invention and biosensors incorporating them. 

DESCRIPTION OF THE FIGURES 

Figure 1 is a schematic diagram of an example of a biosensor of the 
invention. In this embodiment, a self-cleaving DNA is immobilized on a solid matrix 



10 



15 
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that is mounted in a plastic 'spin-column 5 . The self-cleaving DNA remains inactive, 
unless it encounters a specific effector molecule that causes allosteric induction. Test 
sample is added to the porous matrix, allowed to incubate, then the solution is 
collected at the bottom of the tube via centrifugation. Since catalytic activity is a 
5 function of the presence (in concentration) of the effector, the concentration of released 
DNA fragments will report the presence and quantity of effector. 

Figure 2 illustrates a sequence and secondary-structure model for a self- 
cleaving DNA of the invention (SEQ ID NO: 1). The bracket indicates the main 
region of DNA cleavage. 

10 Figure 3 sketches an example of (A) an immobilized DNA biocatalyst of the 

invention and (B) a simple reactor assembly. 

Figure 4 shows a demonstration of catalytic function by immobilized DNA 
enzymes. 5 ' 32 P-labeled RNA substrate was applied to a streptavidin column (AffiniTip 
Strep 20, Genosys Biotechnologies) that was derivatized with 5'-biotinyl DNA 
15 enzyme. The DNA enzyme was immobilized to give an effective concentration of ~ 1 
^iM. Substrate (0.5 jiM was applied to the column in repetitive 20 |aL aliquots, 
allowed to react for 10 min., then recovered for analysis by polyacrylamide gel 
electrophoresis. Fraction of substrate cleaved was plotted as a function of volume 
eluted. 

20 Figure 5 illustrates hammerhead ribozyme constructs described in Example 1 

below. HI (SEQ ID NO: 2) is identical to the ribozyme 'HHJjT that was originally 
characterized by Fedor and Uhlenbeck (12). H2 (SEQ ID NO: 3) carries an additional 
G-C base pair in stem I and is flanked on each end by accessory sequences that are 
designed as short hairpins to reduce the occurance of inactive structures. H3 (SEQ ID 

25 NO: 4) is an integrated hammerhead ribozyme that includes an RNA domain corre- 
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sponding to a truncated version of an ATP- and adenosine-specific aptamer (35). H4 
and H5 are modified versions of H3 that include an aptamer-domain mutation and a 3 
base-pair extension of stem II, respectively. Arrowheads indicate the site of ribozyme- 
mediated cleavage. 

5 Figure 6 shows evidence of ATP- and adenosine-mediated inhibition of a 

hammerhead ribozyme described in Example 1. (A) Hammerhead constructs HI, H2 
and H3 (400 nM) were incubated with trace amounts of (5'- 32 P)-labeled substrate (S) 
in the absence (-) or presence (+) of 1 mM ATP for 30 min. (B) The specificity of 
the effector molecule was examined by incubating H3 and S for 45 min as described 
10 in Example 1 without (-) or with 1 mM of various nucleotides as indicated. Similarly, 
constructs H4 and H5 were examined for activity in the presence of 1 mM ATP. 
Reaction products were separated by a denaturing (8 M urea) 20% polyacrylamide gel 
and visualized by autoradiography. E, S and P identify enzyme, substrate and product 
bands, respectively. 

15 Figure 7 plots kinetic analysis results of the catalytic inhibition of H3 by 

ATP described in Example 1 . (A) Plot of H3 ribozyme activity (400 mM) in the 
presence of 10 jaM (open circles) and 1 mM (filled circles) ATP. Dashed line repre- 
sents the average initial slope obtained in the absence of ATP or in the presence of as 
much as 1 mM dATP. (B) Plot of H3 ribozyme activity (k obs ) in the presence of 

20 various concentrations of dATP (open circles) and ATP (filled circles). Also plotted on 
the y axis are k obs values for HI, H2 and H3 (open squares, filled squares and open 
circles, respectively) with no added effector molecules. 

Figure 8 (A) depicts integrated constructs for allosteric induction by ATP 
(H6, SEQ ID NO: 5 and H7, SEQ ID NO: 6) and allosteric inhibition by theophilline 
25 (H8) described in Example 1 . H7 replaces the central U-A pair with a G e U mismatch 
and is designed to further reduce hammerhead catalysis. H8 is analogous to H3 except 
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that the ATP-aptamer domain is replaced by the theophylline aptamer corresponding to 
*mTCT8-4' that was described by Jenison, et al (21). Arrowhead indicates the site of 
ribozyme-mediated cleavage. (B) Induction of ribozyme catalysis during the course of 
a ribozyme reaction was examined by incubating H6 in the absence (open circles) and 
5 presence (open squares) of 1 mM ATP, and when ATP is added (filled circles), to a 
final concentration of 1 mM during an ongoing ribozyme reaction. Arrow indicates 
the time of ATP addition. 

Figure 9 shows in vitro selection of self-cleaving DNAs described in 
Example 2. In a a y (I) a pool of 5'-biotinylated DNAs is immobilized on a streptavidin 

10 matrix, washed to remove unbound DNAs, then (II) eluted under the desired reaction 
conditions to separate self-cleaving DNAs from those that are inactive. (Ill) Selected 
DNAs are amplified by the polymerase chain reaction (PCR) and (IV) the selection 
round is completed by immobilizing the resulting double-stranded DNAs on new 
matrix followed by removal of the nonbiotinylated strand by chemical denaturation. 

15 (V) The pool is prepared for further analysis by PCR amplification with non-biotiny- 
lated primers. Encircled B indicates 5' biotin. In 6, the construct used for the initial 
round of selection contains a domain of 50 random-sequence nucleotides (N 50 ) flanked 
by 38 and 14 nucleotides of defined sequence. DNAs used in subsequent rounds carry 
an additional 26 nucleotides, as defined by primer 1 (SEQ ID NO: 7; primer 2 is SEQ 

20 ID NO: 8). Precursors that cleave within the overlined region retain sufficient 5' 

primer binding site for amplification and are expected to be favored during selection. 
In c, self-cleavage activity of the initial DNA pool (GO) and the pool isolated after 
seven rounds (G7) of selection. 5' 32 P-labeled precursor DNA (Pre) was incubated in 
the presence (+) or absence (-) of 10 |iM each of Cu 2+ and ascorbate, or in the absence 

25 of Cu 2+ or ascorbate (-Cu 2+ and -asc, respectively) for various times as indicated. M 
is 5' 32 P-labeled primer 3 and Civ identifies cleavage products. 
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Figure 10 shows sequence analysis and catalytic activity of individual G8 
DNAs described in Example 2 (SEQ ID NOs: 9-31). In a, alignment of 34 sequences 
reveal the presence of two major classes of molecules that are characterized by sets of 
common sequences (boxed nucleotides). DNAs that were encountered more than once 
5 are identified by noting the number of occurrences in parentheses. In b 9 self-cleavage 
activity of ~5 nM 5' 32 P-labeled precursor DNA from G8 DNA and from individuals 
CA1, CA2 and CA3 in the absence (-) or presence (+) of Cu 2+ and ascorbate (10 |iM 
each) are shown. 

Figure 1 1 depicts cleavage site analysis of CA3 (lane 2), an optimized 
10 variant (variant 1, Figure 13b) of CA3 (lane 3) and CA1 (lane 4) described in Exam- 
ple 2. DNA size markers (lane 1) are 5' 32 P-labeled DNAs of 10-13 nucleotides as 
indicated. The nucleotide sequence of these markers correspond to the 5' terminal 
constant region of the precursor. 

Figure 12 (a) shows an artificial phylogeny of CA1 (SEQ ID NO: 30) 
15 variants described in Example 2. The numbered sequence is wild-type CA1, and 
nucleotides of variants that differ from this sequence are aligned below. A dash 
indicates a deleted nucleotide, (b) Partial secondary-structure model for a variant of 
CA1 (arrowhead, SEQ ID NO: 32). Numbered nucleotides are derived from the 
region that was randomized in the starting pool. Asterisk indicates the primary 
20 cleavage site and the bar defines the region that undergoes detectable cleavage. Not 
detailed are nucleotides within the 3' primer binding site that are also required for 
catalytic activity. 

Figure 13 shows a Cu 2+ -dependent self-cleaving DNA described in Example 
25 2. (a) Cleavage assay of G8 DNA, CA1 (SEQ ID NO: 30), CA3 (SEQ ID NO: 31) 
and the optimized population of CA3 variants was isolated after mutagenesis followed 
by five additional rounds of selection, (b) Sequence alignment of individual CA3 
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variants that have been optimized for catalytic function with Cu 2+ . The numbered 
sequence is wild-type CA3, and nucleotides of variants that differ from this sequence 
are aligned below. A dash indicates a deleted nucleotide. Arrowheads identify CA3 
variants 1-3 as denoted. Asterisk and bar indicate the major and minor Civ 2 cleavage 
5 sites, respectively. 

Figure 14 shows the sequence and predicted secondary structures of mini- 
mized self-cleaving DNAs described in Example 3. (A) Sequence and secondary 
structure of a synthetic 69-nucleotide self-cleaving DNA that was isolated by in vitro 
selection (SEQ ID NO: 33). Numbers identify nucleotides that correspond to the 50- 

10 nucleotide random-sequence domain that was included in the original DNA pool (note 
that 19 bases of this domain have been deleted). The conserved nucleotides (11-31, 
boxed) are similar to those previously used to define this class of deoxyribozymes 
(Example 2). (B) A 46-nucleotide truncated version of class II DNAs that retains full 
activity (SEQ ID NO: 34). I and II designate stem-loop structures of the 46mer that 

15 are predicted by the structural folding program 'DNA infold' (18, 19), and that were 
confirmed by subsequent mutational analysis (Figure 15). The conserved core of the 
deoxyribozyme spans nucleotides 27-46 and the major site of DNA cleavage is 
designated by the arrowhead. Encircled nucleotides can be removed to create a 
bimolecular complex where nucleotides 1-18 constitute the 'substrate 5 subdomain, and 

20 nucleotides 22-46 constitute the 'catalyst' subdomain. 

Figure 15 shows a confirmation of stems I and II by mutational analysis 
described in Example 3. (A) Trace amounts of 5' 32 P-labeled substrate DNAs (sl-s3) 
were incubated with 5 (iM complementary or non-complementary catalyst DNAs (cl- 
c3) in reaction buffer A containing 10 |iM CuCl 2 at 23°C for 15 min. Reaction 
25 products were separated by denaturing 20% polyacrylamide gel electrophoresis 
(PAGE) and imaged by autoradiography. Bracket identifies the position of the 
substrate cleavage products. (B) Self-cleavage activity of the original 46mer sequence 
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compared to the activity of variant DNAs with base substitutions in stem II. Individual 
46mer variants (100 pM 5' 32 P-labeled precursor DNA) were incubated for the times 
indicated under reaction conditions as described above. Clvl and Clv2 identify 5'- 
cleavage fragments produced upon precursor DNA (Pre) scission at the primary and 
5 secondary sites, respectively. Mutated positions are defined using the numbering 
system given in Figure 14. 

Figure 16 identifies a triplex interaction between substrate and catalyst 
DNAs described in Example 3. A revised structural representation portrays a triple- 
helix interaction (dots) between the four base pairs of stem II and four consecutive 
10 pyrimidine residues near the 5' end of the substrate DNA. c4 and s4 represent 

sequence variants of c3 (SEQ ID NO: 36) and s3 (SEQ ID NO: 35) that retain base 
pairing within stem II, and that use an alternate sequence of base triples. DNA 
cleavage assays were conducted as described in Figure 15,4. Bracket identifies the 
position of the substrate cleavage products. 

15 Figure 17 shows targeted cleavage of DNA substrates using deoxyribozymes 

with engineered duplex and triplex recognition elements. (A) A 101 -nucleotide DNA 
incorporating three different deoxyribozyme cleavage sites was prepared by automated 
chemical synthesis (SEQ ID NO: 37). Each cleavage site consists of an identical 
leader sequence (shaded boxes) followed by a stem I recognition element of unique 

20 sequence. The specific base complementation between the synthetic catalyst DNAs cl, 
c3 and c7 are also depicted. The catalytic core sequences and the leader sequence/stem 
II interactions for each site are identical (inset). Asterisks indicate G-T wobble pairs 
that allow cross reaction between cl and the target for c3. Dots indicate base triple 
interactions. (5) Cleavage of the lOlmer DNA by cl, c3, and c7 was examined by 

25 incubating trace amounts of 5' 32 P-labeled substrate in reaction buffer containing 30 
|liM CuCl 2 at 23°C for 20 min., either in the absence (-) or presence of 5 |iM catalyst 
DNA as indicated. Reaction products were separated by denaturing 10% PAGE and 
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visualized by autoradiography. (Q Similarly, a 100-nucleotide DNA was prepared 
that contained three identical stem I pairing regions (shaded boxes) preceded by eight 
successive pyrimidine nucleotides of unique sequence (SEQ ID NO: 38). Three 
synthetic deoxyribozymes (c9, clO, ell) that carry identical stem I paring elements 

5 (inset) and extended stem II subdomains of unique sequence, were designed to target 
the three cleavage sites exclusively through DNA triplex interactions. (D) Cleavage of 
lOOmer DNA by c9, clO, and ell was established as described in (B) above. Misclea- 
vage is detected for each triplex-guided deoxyribozyme upon extended exposure during 
autoradiography (e.g., ell), indicating that weak-forming triplex interactions allow 

10 some DNA-cleavage activity to occur. 

Figure 18 illustrates the in vitro selection of histidine-dependent deoxyribo- 
zymes described in Example 4. (a) A pool of 4 x 10 13 biotin-modified DNAs was 
immobilized on a streptavidin-derivatized column matrix. Each DNA carries a single 
embedded RNA linkage (rA) and a 40-nucleotide random-sequence domain that is 

15 flanked by regions that are complementary to nucleotides that reside both 5' and 3' of 
the target phosphodiester (pairing elements i and ii; SEQ ID NOs: 39 and 40). These 
pre-engineered substrate-binding interactions are expected to increase the probability of 
isolating active catalysts (7). DNAs that catalyze the cleavage of the RNA linkage 
upon incubation with a solution buffered with histidine were released from the matrix, 

20 were amplified by the polymerase chain reaction (PCR), and the amplification products 
again were immobilized to complete the selection cycle (14-16). (b) Four classes of 
deoxyribozymes were determined by sequence comparison (SEQ ID NOs 41 to 44). 
Variants within each group differed by no more that two mutations from the sequences 
shown. Catalytic assays active (+) when either HEPES or histidine buffers are used, 

25 while class II DNAs not active (-) when histidine is absent. Arrowhead identifies the 
site of cleavage and numbers correspond to the original 40-nucleotide random-se- 
quence domain. 
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Figure 19 shows sequences and secondary structures of variant deoxyribo- 
zymes discussed in Example 4. (a) Individual DNAs isolated after reselection of 
mutagenized pools based on the class II deoxyribozyme (II) (SEQ ID NO: 45) or the 
HD2 deoxyribozyme (HD2 pool, SEQ ID NO: 46). Depicted are the nucleotide 

5 sequences for the mutagenized core of the parent DNAs and the nucleotide changes for 
each variant deoxyribozyme examined after reselection. Deoxyribozymes HD1 (SEQ 
ID NO: 47) and HD2 (SEQ ID NO: 48) were recovered from DNA pools generated 
after five rounds of reselection with 50 or 5 mM histidine, respectively, (b) Each 
deoxyribozyme was reorganized to create a bimolecular complex, whereby separate 

10 substrate molecules are recognized by two regions of base complementation (stems I 
and II) with the enzyme domain. Deoxyribozyme nucleotides are numbered consecu- 
tively from the 5' terminus. 

Figure 20 shows cofactor recognition by a deoxyribozyme described in 
Example 4. (a) Catalytic activity of HD1 with L-histidine, D-histidine, and various 

15 dipeptides that received (+) or did not receive (-) pretreatment with hydrochloric acid. 
HD1 (10 jiM) was incubated in the presence of trace amounts of 5' 32 P-labeled 
substrate oligonucleotide (Figure 19b) and were incubated at 23°C for 2.5 hr with 50 
mM L-histidine, D-histidine, or various dipeptides as indicated. Reaction products 
were analyzed by denaturing (8 M urea) polyacrylamide gel electrophoresis (PAGE) 

20 and imaged by autoradiography. S and P identify substrate and product (5 '-cleavage 
fragment) bands, respectively, (b) Chemical structures of L-histidine and the ana- 
logues used to probe deoxyribozyme cofactor specificity, (c) Representative deoxyri- 
bozyme assays for HD1 (El) catalytic activity with selected amino acids and histidine 
analogues. Reactions and analyses were conducted as described in a. 

25 Figure 21 are graphs showing the involvement of histidine in deoxyribozyme 

function described in Example 4. (a) Concentration-dependent induction of deoxyribo- 
zyme function by histidine. Open and shaded arrowheads indicate the concentration of 
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histidine that was maintained during the selection of HD1 and HD2, respectively, (b) 
Dependence of deoxyribozyme function on pH. Data represented in the main plot was 
produced using 1 mM histidine while data given in the inset was obtained using 5 mM 
histidine. Data depicted with filled, open, and shaded circles was collected using 
5 MES-, Tris-, and CAPS-buffered solutions, respectively. 



DETAILED DESCRIPTION OF THE INVENTION 

Natural ribozymes and artificial ribozymes and deoxyribozymes that have 
been isolated by in vitro selection are not known to operate as allosteric ribozymes. 
This invention is based upon the finding that small-molecule effectors can bind to 

10 ribozyme and deoxyribozyme domains and modulate catalytic rate. As will be dis- 
cussed more fully below, in the practice of the invention, an effector molecule or 
effect binds or affects an allosteric site that is spatially distinct from that of the 
enzyme or reporter domain. Allosteric polynucleotides of the invention can thus 
rapidly interconvert from an "off 1 state to an "on" state, or vice versa, reversibly, on a 

15 time scale that is relevant for their use as biosensors and bioswitches. For example, 
using rational design strategies, a 'hammerhead' self-cleaving ribozyme described 
herein was coupled to different aptamer domains to produce ribozymes whose rates 
can be specifically controlled by adenosine and it's 5'-phosphorylated derivatives. A 
number of other allosteric ribozymes have been created that are sensitive to a variety 

20 of other effectors, including drug compounds, biological metabolites, and toxic metals. 
It is possible to construct, using a mix of in vitro selection and rational design 
strategies, novel biosensors that rely on nucleic acid sensor elements. To achieve this, 
unique RNA or DNA sequences can be appended to ribozymes or deoxyribozymes, 
thereby creating new enzymes having catalytic rates that can be influenced by specific 

25 chemical effectors {e.g., molecules of diagnostic interest), physical signals, and 
combinations thereof. 
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About 50 years ago, it was observed with some polypeptide enzymes that 
catalytic plots of reaction velocity, V, versus substrate concentration [S], displayed 
sigmoidal plots, rather than hyperbolic plots predicted by the simple enzyme + 
substrate model of enzymatic action described by Michaelis-Menten in 1913. In 1965, 

5 Monod, et aL, explained these findings by suggesting that enzymatic reaction rates 
were altered by regulatory domains (3a). In this classical model of "allostery", 
enzymatic activity by "allosteric enzymes" is modulated by reversible binding to com- 
pounds, termed "effectors", at specific sites other than the enzyme's substrate binding 
sites, which, accordingly, are called "allosteric" sites. At constant enzyme and 

10 substrate concentrations, binding of a negative "effector" reduces the reaction rate 

("allosteric inhibition"), and binding of a positive "effector" increases the rate ("alloste- 
ric activation"). Allosteric inhibition may be achieved a number of ways, including 
reducing the binding affinity of the enzyme for its substrate (often reported as increas- 
es in Michaelis-Menton parameter KJ and/or by increasing the time required for each 

15 catalytic turnover (often reported as a decrease in V max ). Conversely, allosteric activa- 
tion may occur either by reduction in or by an increase in V max , or both. 

Decades later it was found that polynucleotides could also catalyze chemical 
reactions, and in 1995, Porta and Lizardi described what they called the first "alloste- 
ric" ribozyme (32a). This was a hammerhead, self-cleaving ribozyme that could be 

20 rendered active by incubating it with a 35-nucleotide antisense DNA oligomer for 

several hours. Notwithstanding the terminology used in the paper, this was not a true 
allosteric effect. Antisense interactions such as that described between the oligonucle- 
otide and the ribozyme are typically comprised of strong base pairing contacts that 
have slow kinetic interchange between bound and unbound states. There was no 

25 allosteric interconversion (from an "off state to an "on" state, or vice versa) disclosed 
upon addition of the 35-mer to an ongoing reaction mixture. Instead, Porta and 
Lizardi described a ribozyme construct which had a folding pathway that could be 
dictated by the 35-mer, but not allosterically switched from active to inactive forms 
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immediately upon addition or depletion of a small effector molecule to or from the 
reaction mixture. Hence, their need for long preincubation and incubation times, and a 
large oligonucleotide that could kinetically and thermodynamically lock the ribozyme 
into an active configuration. 

5 In contrast, in the practice of the invention, purified functional DNA and/or 

polynucleotides that exhibit true allosteric properties that modify a function or 
configuration of the polynucleotide with a chemical effector, a physical signal, or 
combinations thereof, are constructed. The function of polynucleotides of the inven- 
tion is not necessarily controlled by base pairing to an oligonucleotide, but, instead, by 

10 binding of a small molecule effector to an allosteric binding site, or interaction of a 
physical signal with an allosteric site, spatially distinct from the enzyme domain, such 
that the function of the polynucleotide is allosterically modulated. In some embodi- 
ments, the polynucleotide is an enzyme exhibiting allosteric properties that modify the 
rate of catalysis of the enzyme. The invention further provides functional RNA or 

15 DNA polynucleotides having catalytic properties with rates that can be positively 

and/or negatively controlled by a chemical effector, a physical signal, or combinations 
thereof. For example, where enzyme polynucleotides of the invention exhibit a 
reaction rate that is enhanced or inhibited by reversible binding to a chemical effector 
at an allosteric binding site spatially distinct from the substrate binding or self-cleaving 

20 site. In some embodiments, the polynucleotides contain from about 100 or fewer 
bases; others are much larger. 

Allosteric polynucleotides of the invention are comprised of any natural, 
recombinant, or synthetic RNA, DNA and mixtures of RNA and DNA. As used 
herein, the terms M DNA f! and M RNA" specifically include sequences that have RNA 
25 and/or DNA analogues. Analogues include chemically modified bases and unusual 
natural bases. Further encompassed by the invention are polynucleotides modified 
during or after preparation of the domains and constructs using standard means. DNA 
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and/or RNA starting materials for the domains, and constructs and complexes contain- 
ing them, may be isolated from whole organisms, tissues or tissue cultures; construct- 
ed from nucleotides and oligonucleotides using standard means; obtained commercial- 
ly; selected from random and enriched in vitro or in vivo sequence pools; and 
5 combinations thereof. 

Any element, ion, and/or molecule can be used as chemical effectors for 
interaction with the bioreactive allosteric polynucleotides of the invention. It is an 
advantage of the invention that the rational design strategies used to construct the 
polynucleotides (discussed more fully below) can be adapted to a great variety of 

10 effectors. A vast number of ligand-responsive ribozymes with dynamic structural 
characteristics can be generated in a massively parallel fashion (23b). Examples 
include, but are not limited to, organic compounds and mixtures of organic compounds 
and metal ions. Chemical effectors may be amino acids, amino acid derivatives, 
peptides (including peptide hormones), polypeptides, nucleosides, nucleotides, steroids, 

15 sugars or other carbohydrates, pharmaceuticals, and mixtures of any of these. Many 
are small; hence, peptides having 9 or fewer amino acid substitutents and disac- 
charides and trisaccharides are typical polypeptide and carbohydrate effectors. 
Illustrated hereafter are theophylline, ATP and modified ATP; 3-methylxanthine, 
cGMP, cCMP, cAMP, FMN, cobalt, cadmium, nickel, zinc, and manganese have also 

20 been shown to be effectors that modulate the reaction rates of polynucleotides of the 
invention (see, for example, various effectors described in 21a, 23a, 23b, 36a, 39a, 
39b, 39c, and 39e). In many preferred embodiments, small molecule effectors, 
typically having a molecular weight of about 300 or less, are employed, including 
metal ions, amino acids, amino acid derivatives, nucleosides, nucleotides, simple 

25 sugars, and steroids. Effectors can be much larger in other embodiments; larger 
molecule effectors can have molecular weights ranging in the tens or thousands Da, 
and sometimes even larger; protein effectors, for example, can range up to 500,000, 
and sometimes several million, Da. In some embodiments, the chemical effectors are 
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microbial or cellular metabolites or other biological samples. Components found in 
liquid biological samples such blood, serum, urine, semen, tears, and biopsy homogen- 
ates taken from patients for medical or veterinary diagnostic or therapeutic purposes 
are particularly preferred chemical effectors in some embodiments (36a). In industrial 
5 and environmental applications, the effectors are pesticides, herbicides, food toxins, 
product ingredients, reactants, and contaminants, drugs, and the like. Allosteric 
polynucleotides of the invention can be used to detect the presence or absence of 
compounds, as well as their concentration (36a). 

Bioreactive polynucleotides of the invention exhibit allosteric properties that 

10 modify polymer function or configuration with a physical signal or a combination of a 
physical signal and a chemical effector in alternate embodiments. Physical signals 
include, but are not limited to, radiation (23a), temperature changes, movement, 
physical conformational changes in samples, and combinations thereof. Physical 
signals include, but are not limited to, tags, beacons, and the like allosteric reporters 

15 that respond to UV, IR, and/or visible light (23b, 44b). The effects are reversible. 
Chemical effectors binding to allosteric ribozymes and/or deoxyribozymes of the 
invention, for example, can enhance or inhibit the catalytic rate, or do both. It is an 
advantage of the invention that, because the molecules are truely allosteric, any type of 
allosteric interconversion is possible. Hence, a sample of allosteric polynucleotide 

20 enzymes can be fully active, partially active, or fully inactive. In other words, acting 
as a switch, they can be all "on" or all "off, or exhibit any level of activity between 
"on" or "off \ (For a further discussion of switches, see Soukup and Breaker, 39c). 
Morever, because they are truly allosteric, the observable response time to an effector 
molecule or effect is immediate. The kinetics of allosteric polynucleotides are similar 

25 to what is observed with allosteric polypeptides. Illustrated hereafter are polynucleo- 
tides that react in less than 60 minutes, preferably inless than 6 minutes, and most 
preferably, in less than about a minute. (See, for example, Figures 7 and 8.) Most 
preferred allosteric polynucleotides respond to effectors within seconds. 



- 19 - 

Many embodiments employ bioreactive allosteric polynucleotides of the 
invention as biosensors in solution or suspension or attached to a solid support such as 
that illustrated in Figure 1. Alone or as a component of a biosensor, the polynucleo- 
tides are used to detect the presence or absence of a compound or its concentration 
5 and/or a physical signal by contact with the polynucleotide. In a typical practice of 
these methods, a sample is incubated with the polynucleotide or biosensor comprising 
the polynucleotide as a sensing element for a time under conditions sufficient to 
observe a modification or configuration of the polynucleotide caused by the allosteric 
interaction. These are monitored using any method known to those skilled in the art, 

10 such as measurement and/or observation of polynucleotide self-cleavage; binding of a 
radioactive, fluorescent, or chromophoric tag; binding of a monoclonal or fusion 
phage antibody; or change in component concentration, spectrophotometric, or 
electrical properties. It is an advantage of the invention that current biosensor 
technology employing potentiometric electrodes, FETs, various probes, redox media- 

15 tors, and the like can be adapted for use in conjunction with the new polynucleotide 
biosensors of the invention for measurement of changes in polynucleotide function or 
configuration. 

The initial studies described in the Examples that follow have involved the 
creation and characterization of novel RNA- and DNA-cleaving enzymes that function 

20 with specific cofactors, or that can be regulated by specific small-molecule chemical 
effectors, physical signals, or combinations thereof. It is clear that additional mole- 
cules with similar sensor and biocatalytic properties can be created by similar means, 
thereby expanding the applications of such molecules. The creation and characteriza- 
tion of a prototype biosensor for ATP is given herein. One construct (H3) in particu- 

25 lar shows ATP concentration-dependent catalytic activity, indicating that this ribozyme 
could be adapted for use in reporting the concentration of this ligand in test solutions: 
Specifically, H3 RNA actively self-cleaves in concentrations of ATP that are below 1 
micromolar, but is maximally inhibited (170-fold rate reduction) in the presence of 1 
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millimolar ATP (Figure 3b). The catalytic rate of the ribozyme in concentrations of 
ATP that range between these two extremes is reflective of the ATP concentration, and 
can be used to determine unknown concentration values. It is important to note that 
the receptor portion of this allosteric ribozyme is completely artificial (created via in 
5 vitro selection) (35), and could be exchanged for other artificial or natural receptor 
domains that are specific for other ligands. 

New and highly-specific receptors can be made via in vitro selection or 
'SELEX' (4,5) using simple chromatographic and nucleic acid amplification techniques 
(4, and illustrated in the Examples). RNA and DNA 'aptamers' produced in this way 
10 can act as efficient and selective receptors for small organic compounds, metal ions, 
and even large proteins. In a dramatic display of RNA receptor function, a series of 
RNA aptamers for theophilline have been isolated (35) that show ~1 0,000-fold 
discrimination against caffeine, which differs from theophilline by a single methyl 
group. 

15 One can isolate new classes of aptamers that are specific for innumerable 

compounds to create novel biosensors or even controllable therapeutic ribozymes for 
use in medical diagnostics, environmental analysis, etc. In the examples that follow, 
simple design strategies have been used to create conjoined aptamer-ribozyme com- 
plexes who's rates can be controlled by small effector molecules. Preliminary studies 

20 have already shown that theophilline-dependent ribozymes can be created through 

rational design. Theophilline, for example, is an important drug for the treatment of 
asthma and it's therapeutic effect is highly dependent on concentration. A biosensor 
for theophilline concentration would be of significant value. Further examination of 
this allosteric ribozyme and of other model ribozymes will help to lay the biochemical 

25 and structural foundations for the design of additional sensor molecules based on RNA 
and DNA. 
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It is an advantage of the invention that the discovery that DNA can function 
as an enzyme (5) has made practical the engineering of enzymes that are chemically 
more stable than either RNA or proteins. The half-life for the hydrolytic breakdown 
of a DNA phosphoester is -200 million years, making DNA the most stable of the 
5 three major biopolymers. These features of DNA, coupled with the fact that DNA 
also can be made to bind various ligand with great specificity and affinity, make this 
polymer an attractive medium for the creation of new industrial enzymes and as sensor 
elements for diagnostics. Also, modified DNAs can be made that are resistant to 
degradation by natural nucleases, making DNA analogues an attractive format for use 

10 in biological solutions. As illustrated hereafter, it has been found that DNA can be 
made to self cleave in a metal ion-dependent fashion. The creation of these DNAs 
that catalyze their own cleavage in the presence of copper can now be used as a 
sensitive reporter of free copper concentration in solution. Another example given 
below is a polynucleotide reactive to histidine. Further engineering of such catalysts 

15 will yield allosteric DNA enzymes that can be used to detect a wide variety of ligands, 
or that report other reaction conditions such as the concentration of salts, pH, tempera- 
ture, etc. In addition, these DNAs may be conducive to monitoring via amperometric 
H 2 0 2 probes or by spectrophotometry analysis of the redox state of copper. Clearly, 
the diversity of signal read-out for both RNA and DNA sensors can be expanded. 

20 Another feature of the invention is that use of polynucleotides as biosensors 

offer advantages over protein-based enzymes in a number of commercial and industrial 
processes. Problems such as protein stability, supply, substrate specificity and inflexible 
reaction conditions all limit the practical implementation of natural biocatalysts. As 
outlined above, however, DNA can be engineered to operate as a catalyst under defined 

25 reactions conditions. Moreover, catalysts made from DNA are expected to be much more 
stable and can be easily made by automated oligonucleotide synthesis. In addition, DNA 
catalysts are already selected for their ability to function on a solid support and are 
expected to retain their activity when immobilized. 
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The invention further encompasses the use of bioreactive allosteric polynucleo- 
tides attached to a solid support for use in catalytic processes. Immobilizing novel DNA 
enzymes will provide a new form of enzyme-coated surfaces for the efficient catalysis of 
chemical transformations in a continuous- flow reactor under both physiological and non- 
5 physiological conditions. The isolation of new DNA enzymes can be each tailor-made 
to efficiently catalyze specific chemical transformations under user-defined reaction 
conditions. The function of catalytic DNAs to create enzyme-coated surfaces that can be 
used in various catalytic processes is described herein and illustrated in Figure 4. Due to 
the high stability of the DNA phosphodiester bond, such surfaces are expected to remain 
10 active for much longer than similar surfaces that are be coated with protein- or RNA- 
based enzymes. 

A variety of different chromatographic resins and coupling methods can be 
employed to immobilize DNA enzymes. For example, a simple non-covalent method that 
takes advantage of the strong binding affinity of streptavidin for biotin to carry out a 
model experiment is illustrated in Figure 3. In other embodiments, DNA enzymes can 
be coupled to the column supports via covalent links to the matrix, thereby creating a 
longer-lived catalytic support. Various parameters of the system including temperature, 
reaction conditions, substrate and cofactor concentration, and flow rate can be adjusted 
to give optimal product yields. In fact, these parameters can be preset based on the 
kinetic characteristic that are displayed by the immobilized DNA enzyme. However, in 
practice, product formation will be monitored and the chromatographic parameters will 
be adjusted accordingly to optimize the system. 

A prototype system for the large-scale processing of RNA substrates using 
an immobilized DNA enzyme is described herein. Product yields have been deter- 
25 mined by analysis of 32 P-labeled substrate and product molecules by polyacrylamide 
gel electrophoresis of eluant samples. Multiple turn-over of immobilized enzyme 
during tests of the reactive chromatographic surface has been observed (Figure 4). 



15 
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The in vitro selection and engineering of new tailor-made DNA biocatalysts will 
produce catalytic surfaces for practical use and of unprecedented stability and catalytic 
versatility. 

EXAMPLES 

5 The following examples are presented to further illustrate and explain the 

present invention and should not be taken as limiting in any regard. 

Example 1 

As mentioned above, natural ribozymes (8) and ribozymes that have been 
isolated by in vitro selection are not known to operate as allosteric enzymes (6). This 
10 example illustrates allosteric ribozymes. 

Using simple rational design concepts, aptamer domains with hammerhead 
self-cleaving ribozymes (13) were joined in a modular fashion, to create a series of 
catalytic RNAs that are amenable to both positive and negative allosteric control by 

15 small-molecule effectors. Initial efforts were focused on the 40-nucleotide ATP- 
binding aptamer, termed 'ATP-40-1', that was described by Sassanfar and Szostak 
(35). This motif shows a specific affinity for adenosine 5' triphosphate (ATP; K D -10 
|liM) and adenosine, but has no detected affinity for a variety of ATP analogues 
including 2'-deoxyadenosine 5- triphosphate (dATP) or the remaining three natural 

20 ribonucleoside triphosphates. The aptamer also undergoes a significant conformational 
change upon ligand binding, as determined by chemical probing studies. These 
characteristics were exploited to create a conjoined aptamer-ribozyme molecule that 
could be subject to ATP-dependent allosteric control. 
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The initial integrated design, H3, incorporates several key features into an 
otherwise unaltered bimolecular hammerhead ribozyme that is embodied by HI (Figure 
5). Each ribozyme and conjoined aptamer-ribozyme was prepared by in vitro tran- 
scription from a double-stranded DNA template that was produced by the polymerase 
5 chain reaction using the corresponding antisense DNA template and the primers 

5'GAATTCTAATACGACTCACTATAGGCGAAAGCCGGGCGA (SEQ ID NO: 49) 
and 5 GAGCTCTCGCTACCGT (SEQ ID NO: 50). The former primer encodes the 
promoter for T7 RNA polymerase. 50-|il transcription reactions were performed by 
incubating of 30 pmoles template DNA in the presence of 50 mM Tris-HCl (pH 7.5 at 

10 23°C), 15 mM MgCl 2 , 5 mM dithiothreitol, 2 mM spermidine, 2 mM of each NTP, 20 
jaCi (a- 32 P)-UTP and 600 units T7 RNA polymerase for 2 hr at 37°C. RNA products 
were separated by polyacrylamide gel electrophoresis (PAGE), visualized by autora- 
diography and the ribozymes were recovered from excised gel slices by crush-soaking 
in 10 mM Tris-HCl (pH 7.5 at 23°C), 200 mM NaCl and 1 mM EDTA and quantified 

15 by liquid scintillation counting. The RNA substrate was prepared (Keck Biotechnology 
Resource Laboratory, Yale University) by standard solid-phase methods and the 2'- 
TBDMS group was removed by 24-hr treatment with triethylamine trihydrofluoride 
(15 |il per AU 260 crude RNA). Substrate RNA was purified by PAGE, isolated by 
crush-soaking, (5'- 32 P)-labeled with T4 polynucleotide kinase and (y- 32 P)-ATP, and 

20 repurified by PAGE. Even after exhaustive incubation with HI, approximately 45% of 
the RNA remains uncleaved. The kinetic calculations have been adjusted accordingly. 

Superficially, sequences at the 5' and 3' termini were appended to make the 
constructs amenable to amplification by reverse transcription-polymerase chain reaction 
25 methods for future studies. Surveyed independently as H2 (Figure 5), these changes 
causes a 6-fold reduction in k obs compared to HI (rates are summarized in Table 1). In 
addition to the 5 '-and 3'-terinal flanking sequences, H3 includes a modified hammer- 
head stem II that carries the ATP aptamer. The decision to locate the aptamer here 
was made primarily because changes in stem II can have large effects on the catalytic 
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rates of hammerhead ribozymes (28). In the absence of ATP, this alteration causes an 
additional two-fold reduction in rate compared to HI. 

The RNA-cleavage activity of H3 is significantly reduced when incubated 
with 1 mM ATP (Figure 6A). In contrast, ATP has no effect on the cleavage activity 
5 of HI or H2. Moreover, inhibition is observed in the presence of adenosine, but not 
with dATP or the other ribonucleoside triphosphates (Fig 2B). This inhibition is 
highly specific and is consistent with the observed binding specificity of the aptamer 
(35). 

Table 1. Catalytic rates of various ribozyme constructs. Constructs denoted with * 
10 and t, contain either a functional ATP aptamer or a defective ATP aptamer, respec- 
tively. 

kobs(min-l) 

construct stem II none ATP dATP 

0.58 

0.10 

0.054 0.00031 0.053 

0.042 0.061 

0.075 0.13 

0.022 0.12 0.027 

0.0012 0.0098 0.0009 



A AGGCC 

HI J i i I I 

GCC GG 

a AqqcC 
H2 a M M 

G C C G G 

C AAC 

H3* i M I 

GUUG 



H4* 



C A AC . 
MM 
GUUG 



C AAGGCC 

H5 I i i i I I i 

GUUCCGG 



CGUAUGC 

H6* 

GUGUGUG 



CGUGUGC 

H7 I I 

GUGUGUG 
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To investigate the mechanism of inhibition of H3 by ATP, two additional 
integrated constructs (Figure 5) were designed. H4 is identical to H3, but carries a G 
to C point mutation that is expected to eliminate ATP binding by the aptamer domain 
(35). As expected, this mutation eliminates the inhibitory effect of ATP. The 
5 allosteric effect may be due to the proximity of the aptamer and hammerhead domains. 
Specifically, structural models of the hammerhead indicate a parallel orientation for 
stems I and II (32). In the uncomplexed state, the aptamer domain is likely to exist in 
a single or a set of conformational state(s) that allow catalysis to proceed unhindered. 
However, when complexed with ATP, this domain undergoes a conformational change 
10 that presumably causes steric interference between structures that are appended to 

stems I and II. H5 carries an additional three base pairs in helix II, to further separate 
the domains, and is not inhibited by ATP. This is consistent with an allosteric 
inhibition mechanism that involves conformational change and the mutually-exclusive 
formation of aptamer and ribozyme domains. 

15 The inhibitory effect of ATP with H3 has been confirmed and quantitated by 

kinetic analysis. Ribozyme activity assays were conducted with trace amounts of 
substrate and excess ribozyme concentrations that significantly exceed K^. Replicate 
k obs values obtained for HI and H2 at 200, 400 and 800 nM ribozyme concentration 
under identical assay conditions differed by less that two fold, suggesting that for each 

20 construct, k obs values approach V max . Reactions also contained 50 mM Tris-HCl (pH 
7.5 at 23°C) and 20 mM MgCl 2 , and were incubated at 23°C with concentrations of 
effector molecules and incubation times as noted for each experiments. Ribozyme and 
substrate were preincubated separately for -10 min in reaction buffer and also with 
effector molecules when present, and reactions were initiated by combining preincu- 

25 bated mixtures. Assays with H8 were conducted in 50 mM HEPES (pH 7.3 at 23°C), 
500 mM NaCl and 10 mM MgCl 2 . Catalytic rates (k obs ) were obtained by plotting the 
fraction of substrate cleaved versus time and establishing the slope of the curve that 
represents the initial velocity of the reaction by a least-squares fit to the data. Kinetic 
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assays were analyzed by PAGE and were visualized and analyzed on a Molecular 
Dynamics Phosphorimager. When shorter effector-molecule preincubations are used, 
the catalytic burst was more prominent and when encountered, a post-burst slope was 
used in the calculations. Replicate experiments routinely gave k obs values that differed 
5 by less than 50% and the values reported are averages of two or more experiments. 
Equivalent rates were also obtained for duplicate ribozyme and substrate preparations. 

The H3 ribozyme displays different cleavage rates, after a brief burst phase, 
with different concentrations of ATP (Figure 7A), with the curve closely predicting the 
10 K D of the aptamer for its ligand. A plot of k obs versus ATP or dATP concentration 
(Figure 7B) demonstrates that H3 undergoes -170-fold reduction in catalytic rate with 
increasing concentrations of ATP, but is not inhibited by dATP. 

Whether ATP could also be made to function as a positive effector of 
ribozyme function was investigated by designing H6 and subsequently H7 (Figure 8A), 

15 both which were found to display ATP-dependent allosteric induction. H6 is similar to 
H5, except that four Watson/Crick base-pairs in stem II are replaced with less- stable 
G 0 U mismatches. These changes are expected to significantly weaken stem II and 
result in diminished ribozyme activity. It was intended to exploit the fact that the G-C 
pair that begins stem II within the aptamer domain is not paired in the absence of 

20 ATP, but will form a stable pair when ATP is complexed (35), thereby increasing the 
overall stability of the stem and inducing catalytic activity. Indeed, a -5-fold reduc- 
tion in catalytic activity with H6 compared to H5 was found, yet ribozyme function 
could be specifically and fully recovered with ATP. The catalytic rate of H6 is also 
enhanced by ATP when added during the course of the reaction (Figure 8B). 

25 As with allosteric effectors of proteins, there is no true similarity between 

the effector molecule and the substrate of the ribozyme. Substrate and effector occupy 
different binding sites, yet conformational changes upon effector binding result in 
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functional changes in the neighboring catalytic domain. The specificity of allosteric 
control of ribozymes can be exquisite, and in this example the ribozyme activity is 
sensitive to the difference of a single oxygen atom in the effector molecule. 

With similar model studies, a palate of design options and strategic ap- 
5 proaches that can be used to create ribozymes with controlled catalytic activity can be 
built. The principles used here (secondary binding sites, conformational changes, 
steric effects and structural stabilization) as well as others may be generally applicable 
and can be used to design additional allosteric ribozymes, or even allosteric deoxyri- 
bozymes (37). For example, an allosteric hammerhead (H8, Figure 8A) that includes 
10 the theophylline aptamer described by Jenison, et al (21) was designed. This con- 
struct displays a modest 3-fold reduction in ribozyme activity (k obs of 0.006 v. 0.002 
min" 1 ) when theophylline is added to a final concentration of 100 jaM. In addition, 
Sargueil, et al (21) have suggested similar studies with the 'hairpin' self-cleaving 
ribozyme. 

15 Example 2 

The isolation by in vitro selection of two distinct classes of self-cleaving 
DNAs from a pool of random-sequence oligonucleotides are reported in this example. 
Individual catalysts from 'class Y require both Cu 2+ and ascorbate to mediate oxidative 
self-cleavage. Individual catalysts from class II were found to operate with copper as 
20 the sole cofactor. Further optimization of a class II individual by in vitro selection 
yielded new catalytic DNAs that facilitate Cu 2+ -dependent self-cleavage with rate a 
enhancement that exceed 1 million fold relative to the uncatalyzed rate of DNA 
cleavage. 

DNA is more susceptible to scission via depurination/p-elimination or via 
25 oxidative mechanisms than by hydrolysis (27). To begin a comprehensive search for 
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artificial DNA-cleaving DNA enzymes, DNAs that facilitate self-cleavage by a redox- 
dependent mechanism were screened for. Cleavage of DNA by chelates of redox - 
active metals (e.g., Fe 3+ , Cu 2+ ) in the presence of a reducing agent is expected to be a 
more facile alternative to DNA phosphoester hydrolysis due to the reactivity of 
5 hydroxyl radicals that are produced by reduction of H 2 0 2 (z'.e., Fenton reaction). 
Moreover, a variety of natural and artificial 'chemical nucleases' rely on similar 
cleavage mechanisms (38-39). 

Beginning with a pool of ~2 x 10 13 random-sequence DNAs (Figure 13b), 
eight rounds of selection were carried out (5, 10) (see materials and methods section, 

10 below) for DNAs that self-cleave in the presence of CuCl 2 and ascorbate. The DNA 
pool that was isolated after seven rounds (G7 DNA) displays robust self-cleavage 
activity that requires both Cu 2+ and ascorbate (Figure 13c). Trace amounts of non- 
specific DNA cleavage can be detected with Cu 2+ and ascorbate concentrations of 100 
jaM or above, but no cleavage of random-sequence (GO) DNA was detected under the 

15 final selection conditions (10 ^iM of each cofactor). In contrast, incubation of G7 
DNA yields a number of distinct DNA cleavage products, suggesting that the pool 
contains multiple classes of DNAs that promote self-cleavage at unique sites. 

Sequence analysis of individual DNAs from G8 reveals a diverse set of 
catalysts that were divided into two groups (Figure 10a) based on sequence similarities. 

20 Cleavage assays from three representative DNAs (CA1, CA2 and CA3) confirm that 
two distinct classes of catalysts have been isolated (Figure 10b). It was expected that 
the cleavage sites for the selected catalysts would reside exclusively within the first 23 
nucleotides of the original construct (Figure 13b). Cleavage in this region would result 
in release of the molecule from the solid matrix, yet the cleaved molecules would 

25 retain enough of the original primer-binding site to allow amplification by PCR. 

Cleavage elsewhere in a molecule would release a DNA fragment that has lost the 5'- 
terminal primer-binding site, and would be incapable of significant amplification 
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during PCR. Surprisingly, although CA1 promotes DNA cleavage within this expected 
region, CA2 and CA3 each cleave at a primary region (Civ 1) near the 5' terminus as 
expected, and at a distal region (Civ 2) that resides within the domain that was 
randomized in the original DNA pool. The Civ 1/Clv 2 product ratio of CA3 is 
5 approximately 2:1. 

The distribution of cleavage products between the two sites in CA3 is 
expected to result in a significant disadvantage during the selection process. About 
35% of CA3-like molecules cleave within the center of the molecule (and hence are 
probably not amplified), while only about 65% cleave at the expected site and can be 

10 perpetuated in the next round of selection via amplification by PCR. In contrast, 
100% of the catalysts that cleave exclusively in the primer-binding region can be 
amplified, giving individuals from class I an apparent selective advantage. However, 
CA3-like catalysts were found to persist in additional rounds of in vitro selection and 
actually come to dominate the population by generation 13. The success of these 

15 catalysts can be understood, in part, by examining the catalytic rates of CA1 and CA3. 
The cleavage rate (k obs ) of 0.018 min" 1 was obtained for CA1 under the final selection 
conditions, while cleavage at Civ 1 of CA3 occurs with a k obs of 0.14 min* 1 . Despite a 
high frequency of miscleavage, class II catalysts more rapidly cleave at the correct site, 
giving CA3-like catalysts a distinct selective advantage over catalysts from class I. 

20 Cleavage sites for both classes have been further localized by gel-mobility 

analysis of the 5' 32 P-labeled self-cleavage products (Figure 11). CA1 produces a 
major cleavage product with a gel mobility that corresponds to a 9-nucleotide frag- 
ment, and also yields a series of minor products that correspond to DNAs of 3 to 8 
nucleotides. The cleavage site heterogeneity observed for CA1 is consistent with an 

25 oxidative cleavage mechanism that involves a diffusible hydroxyl radical. Typically, 
cleavage of nucleic acids by an oxidative cleaving agent occurs over a range of 
nucleotides, with a primary cleavage site flanked on each side by sites that are cleaved 
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with decreasing frequency. It has been suggested that the frequency of DNA cleavage 
is proportional to the inverse of the distance that separates the target phosphoester 
linkage and the generation site of the hydroxyl radical (18). However, the distribution 
of cleavage products formed by CA1 are indicative of a unique active site that permits 
5 localized DNA cleavage to occur only at nucleotides that immediately flank the 5' side 
of the major cleavage site. 

Similarly, Civ 1 of CA3 consists of a series products that range in mobility 
from 9 to 14 nucleotides, with the major product corresponding to a 12-nucleotide 
DNA (Figure 11). The major product formed upon DNA scission at Civ 2 corre- 

10 sponds to 70 nucleotides, with minor products corresponding to DNAs of 66-69 

nucleotides. The most frequent site of cleavage at Civ 2 is located near position 34 
(G) of the original random-sequence domain. Oxidative cleavage of DNA can proceed 
by a variety of pathways, each that produce distinct cleavage-product termini (22). 
Therefore, conformation of these cleavage sites must now proceed by conducting a 

15 more detailed analysis of the chemical structures of the reaction products. 

To gain insight into the secondary structure of CA1, an artificial phylogeny 
(2) of functional CA1 sequence variants for comparative sequence analysis (47) were 
produced. The 50 nucleotides that corresponds to the original random-sequence 
domain were mutagenized by preparing a synthetic DNA pool such that each wild-type 

20 nucleotide occurs with a probability of 0.85 and each remaining nucleotide occurs with 
a probability of 0.05. The resulting pool was subjected to five additional rounds of 
selection for activity in the presence of 10 |iM each of Cu 2+ and ascorbate. Sequence 
alignment of 39 resulting clones (Figure 12a) reveal two main regions (nucleotides 20- 
28 and 41-50) of strictly-conserved sequence interspersed with regions that tolerate 

25 variation. A total of 25 positions experienced two mutations or less. Other positions 
show sequence covariation, indicating that these nucleotides may make physical 
contact in the active conformation of the deoxyribozyme. For example, A32 and G40 



- 32 - 



frequently mutate to C or T, respectively. This suggests a preference for these bases 
to pair as C-G or A-T. Indeed, this inferred pairing occurs in a region (nucleotides 
28-44) that has considerable base-pairing potential, consistent with the formation of a 
hairpin structure. 

5 Using sequence data and truncation analyses, a partial secondary-structure 

model for CA1 was constructed (Figure 12b). Both the 5'- and 3 '-terminal nucleotides 
show significant base-pairing potential with the substrate domain of the molecule. The 
putative hairpin domain described above (nucleotides 28-44) is flanked by the con- 
served 3' terminus and by a highly-conserved region that is composed mainly of G 

10 residues. It was found that removal of an additional G-rich region that is located in 
the y primer binding site abolishes the catalytic activity of CA1. Extended stretches 
of G residues that form 'G-quartet' structures (46) have been identified in a number of 
other single- stranded DNAs (3,20,26,48). The G-rich sequence in CA1 may also form 
a G-quartet, either independently or with other stretches of G residues that occur 

15 elsewhere in the primary structure of the catalyst. 

CA1 has no detectable activity in the absence of ascorbate, but surprisingly, 
both the G8 population DNA and CA3 display significant cleavage when only Cu 2+ is 
added (Figure 13a). A k obs = 8 x 10" 4 min* 1 for Civ 1 was measured for CA3 in the 
presence of 10 \xM Cu 2+ . In vitro selection was employed to isolate CA3 variants with 

20 enhanced the Cu 2+ -dep endent activity of CA3. CA3 was mutagenized (see above) and 
subjected to five rounds of selection using 10 \\M Cu 2+ as the sole cofactor. Sequence 
alignment of 40 resulting clones (Figure 13b) reveal a single region of highly-con- 
served sequence, spanning nucleotides 15 to 50 of the original random-sequence 
domain. The base identity of 27 nucleotides within this region were found to vary in 

25 three or fewer individuals. The most notable exceptions to this sequence conservation 
are a T deletion between nucleotides 39 and 45, and a T to G mutation that occurs at 
nucleotide 28. In a related selection experiment, active variants of CA3 in which 
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nucleotides 1 through 20 of the original random-sequence domain have been deleted 
were isolated. 

The catalytic activity of the reselected CA3 pool improved by nearly 100- 
fold, with variant DNAs 1, 2 and 3 (Figure 13b) displaying k obs values of 0.052 min" \ 
5 0.033 min" 1 and 0.043 min" 1 , respectively. The uncatalyzed rate of DNA cleavage in 
the presence of Cu 2+ was assessed by incubating 5' 32 P-labeled DNA oligomer (primer 
3) under identical conditions. No Cu 2+ -dependent cleavage of DNA was detected, even 
after a 2-week incubation at 23°C. The overall rate enhancement of the CA3 variants 
was estimated to be considerably greater than 10 6 fold compared to the uncatalyzed 

10 rate. Both CA3 and variant 1 likely proceed via the same DNA cleavage mechanism, 
as evident by their similar catalytic cleavage patterns (Figure 11). A synthetic 87- 
nucleotide version of variant 1 that lacks the 3 '-terminal primer-binding site remains 
active (k ob5 = 0.02 min 1 for Civ 1, 10 \iM Cu 2+ ), while an inhibitory effect is observed 
with 100 \iM Cu 2+ . In addition, the self-cleavage activity of this truncated DNA has a 

15 pH optimum of 7.5, with no specific monovalent cation requirement. Sequential 
deletion of nucleotides from the 5' terminus of this DNA results in a progressive 
reduction in catalytic activity, with a 4-nucleotide deletion resulting in nearly complete 
loss of function. 

The isolation of a variety of self-cleaving DNAs with Cu 2 7ascorbate- 
20 dependence is consistent with an earlier report (23) of site-specific cleavage of a 
single-stranded DNA under similar conditions. These results confirm that DNA is 
indeed capable of forming a variety of structures that promote chemical transforma- 
tions. In addition, the catalytic rates for both classes of self-cleaving DNAs compare 
favorably to those attained by other deoxyribozymes and by natural and artificial 
25 ribozymes. The finding that DNA is also able to perform self-cleavage with Cu 2+ 
alone is unexpected, since the mechanism for the oxidative cleavage of DNA also 
requires a reducing agent such as ascorbate or a thiol compound (38,39). 
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A number of chemical nucleases have been prepared by others and examined 
for their potential as site-specific DNA-cleaving agents. For example, 1,10-phenan- 
throline and similar agents bind DNA, presumably via intercalation, and positions 
copper ions near the ribose-phosphate backbone where formation of a reactive oxygen 
5 derivative favors cleavage of the DNA chain (39). Alternatively, metal-binding 
ligands have been attached to oligonucleotide probes, in order to construct highly- 
specific DNA cleaving agents that recognize DNA by triple-helix formation (26). The 
catalytic DNAs described in this report likely replace the role of chemical nucleases by 
forming their own metal-binding pockets so as to promote region-specific self-cleav- 

10 age. In fact, the addition of 1,10-phenanthroline to a catalytic assay of a synthetic 
class II DNA actually inhibits catalytic function. The optimal Cu 2+ concentration for 
the 87-nucleotide DNA is -10 jaM, with catalytic activity dropping significantly at 
both 1 and 100 \iM Cu 2+ . The inhibitory effect of 1,10-phenanthroline might be due 
to the reduction in concentration of free Cu 2+ upon formation of Cu 2+ -phenanthroline 

15 complexes. 



While not wishing to be bound to any theory, several different mechanisms 
for the oxidative cleavage of class II DNAs seem possible. For example, the class II 
DNAs may simply scavenge for trace amounts of copper and reducing agents that are 
present in the reaction buffer. Alternatively, these DNA molecules might make use of 

20 an internal chemical moiety as the initial electron donor. In each example, the 

catalytic DNAs could still cleave by an oxidative mechanism, but would at least appear 
to gain independence from an external source of reducing agent. The importance of 
H 2 0 2 in oxidative processes can be examined with catalase, an enzyme that efficiently 
promotes the dismutation of H 2 0 2 molecules to yield water and molecular oxygen. The 

25 catalytic activity of a representative DNA from class II is completely inhibited upon 
the addition of catalase, consistent with the notion that H 2 0 2 is a necessary intermedi- 
ate in an oxidative pathway for DNA cleavage. The catalytic rate of CA3 variants is 
greatly increased when incubated in the presence of added H 2 0 2 . For example, the 87- 



- 35 - 



nucleotide DNA can be made to cleave quantitatively at Civ 1 (k obs = 1.5 min ') in the 
presence of 10 \xM Cu 2+ and 35 mM H 2 0 2 . 

It has not been determined whether trace amounts of H 2 0 2 in water are used 
by the catalysts, or if the DNA can produce H 2 0 2 in the absence of a reducing agent. 
5 It was found that preincubation of separate solutions of catalytic DNA in reaction 
buffer (minus Cu 2+ ) and of aqueous Cu 2+ , followed by thermal denaturation of the 
catalase, results in full self-cleavage activity upon mixing of the two solutions. We 
also find that self-cleavage of the 87-nucleotide variant reaches a combined maximum 
(Civ 1 + Civ 2) of ~ 70%, regardless of the concentration of catalytic DNA present in 

10 the reaction. Similarly, preincubation of a reaction mixture with excess unlabeled 

catalyst (1 jaM) followed by the addition of a trace amount of identical 5' 32 P-labeled 
catalysts produces normal yields of labeled-DNA cleavage products. Finally, addition 
of fresh reaction buffer to a previously-incubated reaction mixture does not promote 
further DNA cleavage, as might be expected if limiting amounts of reducing agent 

15 were responsible for activity. 

Certain constructs of the self-splicing ribozyme of Tetrahymena have been 
shown to catalyze the cleavage of DNA via a transesterification mechanism (19,33), 
and the ribozyme from RNase P has been found to cleave DNA by hydrolysis (31). 
Such ribozymes might also be made to serve as therapeutic DNA-cleaving agents, 

20 analogous to the function of RNA-cleaving 'catalytic antisense' ribozymes (9). The 
secondary-structure model of CA1 (Figure 12b) includes stretches of predicted base 
pairing both 5' and 3' to the primary cleavage site, suggesting that 'substrate' and 
'enzyme' domains can be separated. Likewise, preliminary analysis of class II mole- 
cules reveals similar base complementation. It is expected that both class I and class 

25 II DNAs can be engineered to create catalytic DNAs that specifically cleave DNA 
substrates with multiple turn-over kinetics. 
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In summary, two distinct classes of DNAs that promote their own cleavage 
have been isolated. One class requires copper and catalyzes the oxidative cleavage of 
DNA with a rate in excess of 1 million fold. Extensive regions of both classes of self- 
cleaving DNAs are important for the formation of catalytic structures, as implicated by 
5 sequence conservation found with selected individuals. These results support the view 
that DNA, despite the absence of ribose 2'-hydroxyl groups, has considerable potential 
to adopt higher-ordered structures with functions that are similar to ribozymes. 

Materials and Methods 

Oligonucleotides 

10 All synthetic DNAs were prepared by automated chemical synthesis (Keck 

Biotechnology Resource Laboratory, Yale University). The starting pool is composed 
of DNAs that carry a 5 '-terminal biotin moiety and a central domain of 50 random- 
sequence nucleotides. Primer 3 is an analogue of primer 1 (Figure 1 3b) that contains a 
3 '-terminal ribonucleoside. Primer 4 is the nonbiotinylated version of primer 2 (Figure 

15 13b). Primer 5 is the 5'-biotinylated form of primer 1. 

In vitro selection 

A total of 40 pmoles of pool DNA in 40 |il buffer A (50 mM HEPES, pH 
7.0 at 23°C, 0.5 M NaCl, 0.5 M KC1) was loaded on two streptavidin-matrix columns 
(Affinitip Strep20, Genosys Biotechnologies) and incubated for -5 min. Unbound 

20 DNAs were subsequently removed from each column by pre-elution with 500 jal of 
buffer A, then by 500 (Jil 0.2 N NaOH, and the resulting matrix-bound DNAs were 
equilibrated with 500 ^il buffer A. Catalytic DNAs were eluted with three successive 
20-|il aliquots of buffer B (buffer A, 100 jiM CuCl 2 , 100 |iM ascorbate) for rounds 1- 
3, or buffer C (buffer A, 10 jaM CuCl 2 , 10 |iM ascorbate) for rounds 4-8. Eluate from 

25 each column was combined with 120 jil 4 mM EDTA and 40 pmoles each of primers 
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1 and 2. Selected DNAs and added primers were recovered by precipitation with 
ethanol and amplified by PCR a 200 |il reaction containing 0.05 U jal-1 Taq polymer- 
ase, 50 mM KC1, 1.5 mM MgCl 2 , 10 mM Tris-HCl (pH 8.3 at 23°C), 0.01 % gelatin, 
and 0.2 mM each dNTP for 25 cycles of 10 sec at 92°C, 10 sec at 50 °C and 30 sec at 
5 72°C. The 5 '-terminal region of each cleaved DNA, including the biotin moiety, was 
reintroduced at this stage. Subsequent rounds were performed by immobilizing 20 
pmoles of pool DNA on a single streptavidin column and selected DNAs were 
amplified in a 100 \x\ reaction for 10 to 20 temperature cycles. Steps II-IV (Figure 
13) were repeated until the population displayed the desired catalytic activity, at which 
10 time the pool was PCR amplified with primers 1 and 3, cloned (Original TA Cloning 
Kit, Invitrogen) and sequenced (Sequenase 2.0 DNA Sequencing Kit, U. S. Biochemi- 
cals). Reselections with CA1 and CA3 were initiated with 20 pmoles synthetic DNA. 
This is expected to offer near comprehensive representation of all sequence variants 
with seven or fewer mutations relative to wild type. 

15 Catalytic assays 

5'- 32 P-labeled precursor DNA was prepared by PCR-amplifying double- 
stranded DNA populations or plasmid DNA using 5'- 32 P-labeled primer 4 and either 
primer 5 or primer 3. The antisense strand is removed either by binding the bio- 
tinylated strand to a streptavidin matrix (primer 5) or by alkaline cleavage of the RNA 

20 phosphodiester-containing strand, followed by PAGE purification (primer 3). DNA 

self-cleavage assays (~5 nM 5' 32 P-labeled precursor DNA) were conducted at 23°C in 
buffer A, with cofactors added as detailed for each experiment. For both in vitro 
selection and for assays, reaction buffers that contained ascorbate were prepared just 
prior to use. Self-cleavage assays conducted with catalase (bovine liver, Sigma) 

25 contained 50 mM HEPES (pH 7.0 at 23°C), 50 mM NaCl, 10 \xU CuCl 2 , and 0.5 U/jil 
catalase, and were incubated at room temperature for 20 min. Catalase activity was 
destroyed by heating at 90°C for 5 min. Products were separated by denaturing (8 M 
urea) polyacrylamide gel electrophoresis (PAGE) using a 10% gel and visualized by 
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autoradiography or visualized and quantitated by Phosphorlmager (Molecular Dynam- 
ics). 

Cleavage product analysis 

Primary cleavage sites for CA1 and CA3 were identified by incubating 5' 
5 32 P-labeled precursor DNA in buffer C and assessing the gel mobility of the 5'- 

terminal cleavage fragments by analysis using a denaturing 20% PAGE as compared to 
a series of 5' 32 P-labeled synthetic DNAs that correspond in sequence to the 5' 
terminus of the precursor DNAs. Products resulting from sission at Civ 2 were 
analyzed by denaturing 6% PAGE. 

10 Kinetic analysis 

Catalytic rates were obtained by plotting the fraction of precursor DNA 
cleaved versus time and establishing the slope of the curve that represents the initial 
velocity of the reaction as determined by a least-squares fit to the data. Kinetic assays 
were conducted in buffer C or in buffer A plus 10 |iM CuCl 2 as indicated for each 

15 experiment. Rates obtained from replicate experiments differed by less than two fold 
and the values reported are averages of at least two analyses. 



Example 3 

This example describes a DNA structure that can cleave single-stranded 
DNA substrates in the presence of ionic copper. This deoxyribozyme can self-cleave, 
20 or it can operate as a bimolecular complex that simultaneously makes use of duplex 
and triplex interactions to bind and cleave separate DNA substrates. DNA strand 
scission proceeds with a k obs of 0.2 min \ a rate that is ~10 ,2 -fold faster than the 
uncatalyzed rate of DNA phosphoester hydrolysis. The duplex and triplex recognition 
domains can be altered, making possible the targeted cleavage of single-stranded 
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DNAs with different nucleotide sequences. Several small synthetic DNAs were made 
to function as simple 'restriction enzymes' for the site-specific cleavage of single- 
stranded DNA. 

A Minimal Cu 2+ -Dependent Self-cleaving DNA. In Example 2, a variety of 
5 self-cleaving DNAs were isolated by in vitro selection from a pool of random-se- 
quence DNAs. Most individual DNAs that were isolated after eight rounds (G8) of 
selection conformed to two distinct classes, based on similarities of nucleotide 
sequence and DNA cleavage patterns. Although individual DNAs from both class I 
and class II require Cu 2+ and ascorbate for full activity, the G8 DNA population 

10 displays weak self-cleavage activity in the presence of Cu 2+ alone. A representative 
class II DNA termed CA3 was further optimized for ascorbate-independent activity by 
applying in vitro selection to a DNA pool that was composed of mutagenized CA3 
individuals. The sequence data from this artificial phylogeny of DNAs indicates that 
as many as 27 nucleotides, most of them located near the 3' terminus of the molecule, 

15 are important for self-cleavage activity. 

Beginning with the original G7 DNA population, an additional six rounds of 
in vitro selection was carried out for DNAs that self-cleave in the presence of 10 |iM 
Cu 2+ , without added reducing agent. Analysis of the G13 population of DNAs 
revealed robust self-cleavage activity, demonstrating that catalytic DNAs can promote 
20 efficient cleavage of DNA using only a divalent metal cofactor. The G13 population 
displays the same cleavage pattern that was observed with individual class II DNAs, 
indicating that class II-like DNAs dominate the final DNA pool. 

A total of 27 individual DNAs from G13 were sequenced and, without 
exception, each carried a 21 -nucleotide sequence domain that largely conformed to the 
25 consensus sequence that was used previously to define class II self-cleaving DNAs. 
Although individuals that have a strictly conserved core (spanning nucleotides 1 1 to 
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31, Figure \4A) dominate the G13 pool, two common variations from this consensus 
sequence include a C to T mutation at position 17 (6 of 28 individuals) or the presence 
of six successive Ts instead of five in the region spanning nucleotides 21 to 25 (4 of 
27 individuals). However, significant differences in nucleotide sequence were found to 

5 occur outside this conserved domain, indicating that large portions of the class II 
deoxyribozymes isolated may not be necessary for catalytic activity. Indeed, three 
individual DNAs were found to have undergone deletions of 16, 19, and 20 nucleo- 
tides within the 50-nucleotide domain that was randomized in the original starting 
pool. The predicted secondary structure for the 1 9-nucleotide deletion mutant (69mer 

10 DNA, Figure 144), obtained by the Zucker 'DNA infold* program (33,50; the DNA 
mfold server can be accessed on the internet at www.ibc.wustl.edu/~zuker/dna/forml. 
cgi.), indicates the presence of three base-paired regions; two involve pairing between 
the original random-sequence domain and the 'substrate' domain, and one that involves 
putative base-pairing of nucleotides that lie within the conserved-sequence region. A 

15 synthetic DNA corresponding to the 69-mer depicted in Figure 14,4 undergoes Cu 2+ - 
dependent self-cleavage at two locations with a combined catalytic rate of approxi- 
mately 0.3 min" 1 under the conditions used for in vitro selection (see Materials and 
Methods below for additional discussion on catalytic rates). 

Whether the two pairing regions of the 69-mer that lie within the variable- 
20 sequence region could be replaced by a smaller stem-loop structure was tested by 
synthesizing a 46-mer DNA, in which 26 nucleotides of this imperfect hairpin were 
replaced by the trinucleotide loop GAA (Figure 145). As expected, the truncated 
4 46mer' DNA retains full catalytic activity, thereby confirming that the deleted 
nucleotides are not essential for deoxyribozyme function. This 46-nucleotide deoxyri- 
25 bozyme is predicted to adopt a pistol-like secondary structure (Figure 145) composed 
of two base-paired structural elements (stems I and II) flanked by regions of single- 
stranded DNA. The primary site of DNA cleavage is located at position 14 which 
resides within one of the putative stem structures of the 46mer. The catalyst also 
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promotes DNA cleavage within a region located apart from the main cleavage site 
(Example 2), as might be expected for a deoxyribozyme that makes use of an oxida- 
tive cleavage mechanism (22). 

Bimolecular Deoxyribozyme Complexes: Substrate Recognition by Duplex 
5 and Triplex Formation. Separate 'substrate' and 'catalyst' DNAs can be created from 
the 46mer by eliminating the connecting loop of stem I (Figure 145). Active bimolec- 
ular complexes then can be reconstituted by combining independently prepared 
substrate and catalyst DNAs. Both the unimolecular 46mer and the bimolecular 
complexes examined cleave with identical rates, promoting primary-site cleavage with 

10 a k ohs of approximately 0.2 min* 1 . The importance of stem I was confirmed (Figure 
\5A) by synthesizing different catalyst DNAs (cl, c2 and c3) and assessing their 
ability to cleave different substrate molecules (si, s2 and s3). For example, cl displays 
activity with its corresponding substrate (si), but not when the non-complementary 
substrate DNAs s2 or s3 are substituted. Likewise, c2 and c3 only cleave their 

15 corresponding substrate DNAs s2 and s3, respectively. Extending stem I to create a 
more stable interaction was also found to confer greater binding affinity between 
substrate and catalyst oligonucleotides. These data indicate that base pairing interac- 
tions that constitute stem I are an essential determinant for catalyst/substrate recogni- 
tion. 

20 Stem II was examined by a similar approach using mutant versions of the 

46mer self-cleaving DNA. A series of variant deoxyribozymes with one or two 
mutations included in the putative stem II structure were synthesized and assayed for 
catalytic activity (Figure 155). Disruption of the original C35-G43 base pair in stem 
II, either by mutation of C to G at position 35 or mutation of G to C at position 43, 

25 results in a substantial loss of activity. Cleavage activity is partially restored when 
these mutations are combined in the same molecule to produce a G35-C43 base pair. 
These results are consistent with the stem-loop structure modeled in Figure 14. Addi- 
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tional support for the presence of stem II was found upon sequence analysis of the 
deoxyribozymes that are present in the original in vrtro-selected pool of DNAs. A 
single self-cleaving DNA was found with a core sequence that differs significantly 
from that of the most frequently represented deoxyribozyme. Nucleotides 38-40 of the 
5 more common 46mer sequence are replaced in the variant deoxyribozyme with the 
nucleotides 5-CTGGGG. This alternative sequence extends stem II by a single C-G 
base pair, consistent with the formation of the predicted stem-loop structure. 

Although the existence of stem II is supported by the data derived from 
mutational analysis, the fact that total restoration of deoxyribozyme activity was not 

10 achieved with restoration of base complementation indicates that the identities of the 
base pairs in this structural element are important for maximal catalytic function. 
Moreover, it was found that mutation or deletion of nucleotides 1 -7 of the 46mer 
result in a dramatic loss of DNA cleavage activity. It was recognized that nucleotides 
4-7 within this essential region of the substrate form a polypyrimidine tract that is 

15 complementary to the paired sequence of stem II for the formation of a YR*Y DNA 
triple helix (14). 

To examine the possibility of triplex formation in the active structure of the 
deoxyribozyme, we modified both the base pairing sequence of stem II (c4) and the 
sequence of the polypyrimidine tract of the substrate (s4) to alter the specificity, yet 

20 retain the potential for forming four contiguous base triples (Figure 16). The c4 

variant DNA cleaves its corresponding s4 DNA substrate, but shows no activity with a 
substrate that carries the original polypyrimidine sequence. It was found that even 
single mutations within stem II (e.g., Figure 155) or single mutations within the poly- 
pyrimidine tract cause significant reductions in catalytic activity. However, the 

25 introduction of six mutations in a manner that is consistent with triplex formation 
results in a variant (c4/s4) complex that displays full DNA cleavage activity. This is 
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the first example of a catalytic polynucleotide, natural or artificial, that makes use of 
an extended triple helix for the formation of its active structure (43). 

Targeted Cleavage of DNA 'Restriction Sites ' with Deoxyribozymes. The 
results described above demonstrate that class II deoxyribozymes identify substrate 
5 DNAs by simultaneously utilizing two distinct recognition domains that are formed 
separately by stems I and II. These structures might be further exploited as recogni- 
tion elements to engineer deoxyribozymes that selectively cleave DNAs at different 
target sites. To demonstrate this capability, a 101 -nucleotide DNA that carries three 
identical leader sequences, each followed by different stem I recognition sequences 

10 was synthesized (Figure 17^4). Three catalyst DNAs (cl, c3 and c7) each were 

designed to be uniquely complementary to one of the three target sites. When incubat- 
ed separately with lOlmer substrate, DNAs c3 and c7 cleave exclusively at their 
corresponding target sites, while cl cleaves at its intended site and also to a lesser 
extent at the c3 cleavage site (Figure 175). The cross reactivity observed with cl can 

15 be explained by examining the base-pairing potential of stem I. Of the six nucleotides 
in the cl recognition sequence, four can form standard base pairs, while the remaining 
two form G-T wobble pairs. The contribution of both duplex and triplex recognition 
elements presumably allows for detectable cleavage activity at this secondary location. 

The triplex interaction that is defined by the base-pairing sequence of stem 
20 II can also be exploited to target specific DNA substrates. We designed three new 
catalyst DNAs (c9, clO and ell) that carry identical stem I pairing subdomains, but 
that have expanded and unique stem II subdomains (Figure 17C). When incubated 
separately with a 1 00-nucleotide DNA that carries three uniquely complementary 
polypyrimidine tracts, each catalyst DNA cleaves its corresponding target site with a 
25 rate that corresponds well with that found for the original self-cleaving DNA. In this 
example, substrate selectivity is determined almost entirely by triplex formation, 
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despite the presence of identical and extensive base complementation (stem I) between 
catalyst and substrate molecules. 

Although DNA cleavage catalyzed by the deoxyribozyme is focused within 
the substrate domain, substantial (-30%) cleavage occurs within the conserved core of 

5 the catalyst strand. This collateral damage causes inactivation of the deoxyribozyme 
and, as a result, super-stoichiometric amounts of catalyst DNA are needed to assure 
quantitative cleavage of DNA substrate. Cleavage of the substrate subdomain proceeds 
more rapidly than does cleavage within the catalytic core. In the presence of excess cl, 
si is cleaved at a rate of approximately 0.2 min' 1 (reaction buffer containing 30 jaM 

10 CuCl 2 ), reaching a plateau of -80% cleaved after 20 min. In contrast, cleavage of cl 
in the presence of excess si proceeds more than 2-fold slower, consistent with our 
earlier report that the ratio of self-cleavage localized in the substrate domain to self- 
cleavage in the catalytic core gives a ratio of -2:1. It was established that, barring 
inactivation by miscleavage, the catalyst strand can undergo multiple turnover. 

15 Cleaving Double-stranded DNA by Thermocy cling. Class II catalyst DNAs 

are not able to cleave target DNAs when they reside within a duplex. The catalyst 
DNA, with its short recognition sequence, presumably cannot displace the longer and 
more tightly-bound complementary strand of the target in order to gain access to the 
cleavage site. It was found that an effective means for specific cleavage of one strand 

20 of an extended DNA duplex makes use of repetitive cycles of thermal denaturation and 
reannealing. For example, c3 remains inactive against a double-stranded DNA target 
in the absence of thermal cycling, but efficiently cleaves the same DNA substrate upon 
repeated heating and cooling cycles. Cleavage of the radiolabeled target is quantitative 
after 6 thermal cycles. DNA cleavage by class II DNAs occurs within the base-pairing 

25 region corresponding to stem I, presumably when this region is in double-helical form. 
This, coupled with the observation of substrate recognition by triplex formation, 
suggests that different DNA enzymes might be engineered to cleave duplex DNA 
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substrates without the need for thermal denaturation. Such deoxyribozyme activity 
would be similar to that performed by a number of triplex-forming oligonucleotides 
that have been engineered to bind and cleave duplex DNA using a chemically-tethered 
metal complex such as Fe-EDTA (24-27). 

5 Conclusions. In its unimolecular arrangement, the class II deoxyribozyme 

could be used to confer the capacity for self-destruction to an otherwise stable DNA 
construct. In its bimolecular form, the deoxyribozyme can act as an artificial restriction 
enzyme for single-stranded DNA, whereas protein-based nucleases that cleave non- 
duplex DNA do not demonstrate significant sequence specificity. It is likely that 

10 Ymaximal discrimination by class II catalysts between closely related target sequences 
can be achieved through careful design of the duplex and triplex recognition domains. 
This is expected to eliminate the cross reactivity that was observed here. Although the 
role of most nucleotides within the substrate domain are involved in substrate recogni- 
tion, the importance of each nucleotide within the leader sequence has yet to be fully 

15 delineated. However, guided by the basic rules of duplex and triplex formation, one 
w3can now engineer highly-specific deoxyribozymes that can catalyze the cleavage of 
single-stranded DNA at defined locations along a polynucleotide chain. 



Materials and Methods 

Oligonucleotides 

20 Synthetic DNAs were prepared by automated chemical synthesis (Keck 

Biotechnology Resource Laboratory, Yale University), and were purified by denaturing 
(8M urea) polyacrylamide gel electrophoresis (PAGE) prior to use. Double-stranded 
lOlmer DNA was prepared by the polymerase chain reaction (PCR) as described in 
Example 2 using the primer DNAs 5' 32 P-GTCGACCTGCGAGCTCGA, (SEQ ID NO: 
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51) 5 'GTAGATCGTAAAGCTTCG (SEQ ID NO: 52) and the lOlmer DNA oligomer 
(Figure \7A) as template. 

In yitro selection 

Optimization of class II self-cleaving DNAs was achieved by in vitro 
5 selection essentially as described In Example 2 using a reaction mixture for DNA 
cleavage composed of 50 mM HEPES (pH 7.0 at 23°C), 0.5 M NaCl, 0.5 M KC1 
(buffer A), and that included 10 jaM CuCl 2 . The selection process was initiated with 
20 pmoles G7 PCR DNA in which the 5' terminus of each catalyst strand carried a 
biotin moiety, thereby allowing DNA from this and subsequent generations to be 
10 immobilized on a streptavidin-derivatized chromatographic matrix. Reaction time was 
15 min. for immobilized DNA from G8-G10 and 12, 7 and 5 min. for the G11-G13 
DNA populations, respectively. Individual self-cleaving DNAs from G13 were 
analyzed by cloning (Original TA Cloning Kit, Invitrogen) and sequencing (Sequenase 
2.0 DNA Sequencing Kit, U.S. Biochemicals). 

15 DNA cleavage assays 

To assess the DNA cleavage activity of self-cleaving molecules, radiolabeled 
precursor DNA was prepared by enzymatically tagging the 5 ' terminus of synthetic 
DNAs in a reaction containing 25 mM CHES (pH 9.0 at 23°C), 5 mM MgCl 2 , 3 mM 
DTT, 1 nM DNA, 1.2 jaM (y- 32 P)-ATP (-130 >xCi), and 1 U/^L T4 polynucleotide 

20 kinase, which was incubated at 37°C for 1 hr. The resulting 5' 32 P-labeled DNA was 
isolated by denaturing PAGE and recovered from the gel matrix by crush-soaking in 
10 mM Tris-HCl (pH 7.5 at 23°C), 0.2 M NaCl, and 1 mM EDTA. The recovered 
DNA was concentrated by precipitation with ethanol and resuspended in deionized 
water (Milli-Q, Millipore). Self-cleavage assays using trace amounts of radiolabeled 

25 precursor DNA (-100 pM) were conducted at 23°C in buffer A containing CuCl 2 as 
indicated for each experiment. Examinations of the DNA cleavage activity of bimo- 
lecular complexes were conducted under similar conditions using trace amounts of of 
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5' 32 P-labeled 'substrate' DNA. Cleavage products were separated by denaturing 
PAGE, imaged by autoradiography or by Phosphorlmager (Molecular Dynamics) and 
product yields were determined by quantitation (ImageQuant) of the corresponding 
precursor and product bands. 

5 Kinetic analyses 

Catalytic rates were estimated by plotting the fraction of precursor or 
substrate DNA cleaved versus time and establishing the slope of the curve that 
represents the initial velocity of the reaction as determined by a least-squares fit to the 
data. Upon close examination, DNA cleavage in both the substrate and enzyme 
10 domains displayed a brief lag phase that complicates the determination of the initial 
cleavage rate. In order to avoid the lag phase, the initial slope was calculated only 
using data collected after the reaction had proceeded for 1 min. Rates obtained from 
replicate experiments differed by less than 50% and the values reported are averages of 
at least three analyses. 



15 Example 4 

The in vitro selection of a catalytic DNA that uses histidine as the active 
component for an RNA cleavage reaction is described in this example. An optimized 
deoxyribozyme only binds to L-histidine or to several closely-related analogues and 
subsequently catalyzes RNA phosphoester cleavage with a rate enhancement of -10- 
20 million fold over the uncatalyzed rate. While not wishing to be bound to any theory, 
the DNA-histidine complex apparently performs a reaction that is analogous to the first 
step of the catalytic mechanism of RNase A, in which the imidazole group of histidine 
acts as a general base catalyst. 
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The class of deoxyribozymes that catalyze the cleavage of an RNA phospho- 
ester bond using the amino acid histidine as a cofactor described herein is depicted in 
Figure 1 8a. To assure that metal-dependent deoxyribozymes were not recovered from 
the random-sequence pool of DNAs, the divalent metal-chelating agent ethylenedi- 

5 mainetetraacetic acid (EDTA) was included in a reaction mixture that was buffered 
with 50 mM histidine (pH 7.5). After 11 rounds of selective amplification, the DNA 
pool displayed RNA phosphoester-cleaving activity, both under in vitro selection 
conditions, and in a reaction buffer containing HEPES (50 mM, pH 7.5) in place of 
histidine. Individual molecules cloned from the final DNA pool were grouped into 

10 one of four sequence classes (Figure 19b), and representative clones were tested for 
catalytic activity. Only class II DNAs demonstrate complete dependence on histidine 
while the remaining classes appear to operate independently of any metal ion or small 
organic cofactor. 

The catalytic rate for the original class II deoxyribozyme was ~ 1000-fold 
15 slower (k ohs = 1.5 x 10" 3 min" 1 ) than most natural self-cleaving ribozymes (44). As a 
result, further optimization of catalytic activity was sought in order to provide an 
artificial phylogeny of variant catalysts for comparative sequence analysis. A new 
DNA pool was prepared based on the sequence of class II deoxyribozymes, such that 
the 39 nucleotides corresponding to the original random-sequence domain were 
20 mutagenized with a degeneracy of 0.21 (6). Beginning with a mutagenized pool that 
sampled all possible variant DNAs with seven or fewer mutations relative to the 
original class II sequence, parallel reselection was conducted using reaction solutions 
buffered with either 50 mM histidine, or with 5 mM histidine and 50 mM HEPES. 
Individual DNAs isolated from the populations resulting from five rounds of reselec- 
25 tion are more active than the original class II deoxyribozyme, and show specific 
patterns of conserved sequences and mutation acquisition (Figure 19a). 
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It was speculated that engineered pairing element i included in the original 
DNA construct (Figure 18a) was being utilized by class II deoxyribozymes. In 
contrast, it was recognized that a conserved-sequence domain near the 3 ' end of the 
core (Figure 19a, nucleotides 32-36) was identical to pairing element ii . Considering 
5 these observations, individual deoxyribozymes HD1 and HD2 were designed to operate 
as separate substrate and enzyme domains (Figure 19b). Specificity for the substrate 
oligonucleotide is defined by the Watson/Crick base complementation between the 
substrate and the two pairing arms of the enzyme domain. Class II deoxyribozymes 
have an absolute requirement for histidine as show by the activity of the bimolecular 
10 HD1 construct to 'caged 5 histidine delivered in the form of dipeptides, and to free 
amino acids that were liberated from each dipeptide by acid hydrolysis (Figure 16a). 
In addition, HD1 accepts L-, but not D-histidine as a cofactor. However, samples of 
D-histidine become active upon treatment with HC1 in accordance with the accelerated 
rate of interconversion between the two isomeric forms in acidic conditions (11). 

15 A larger panel of histidine analogues were examined (Fogire 16b) in order to 

more carefully examine the chemical groups of histidine that are important for 
catalytic activity and to rule out the possibility that catalysis might be due to a 
contamination of a metal ion cofactor. HD1 discriminates against a variety of histidine 
analogues, but shows full activity with the methyl ester of L-histidine (Figure 16c). 

20 Both the 1 -methyl- and 3-methyl-L-histidine analogues do not support HD1 activity, 
indicating that the imidazole ring of histidine is important for deoxyribozyme function. 
As expected, HD2 has a similar pattern of cofactor discrimination (Table 2). Both 
catalysts show stereospecific recognition of histidine, and make use of interactions with 
the oc-amino group, with both carboxyl oxygens, and with the imidazole group in 

25 order to attain maximize cofactor binding. Although a number of analogues cannot 
support deoxyribozyme activity, no compounds function as competitive inhibitors, 
indicating that their inactivity is due to the failure to bind the deoxyribozyme. 
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Table 2. Relative k obs values for HD2 in the presence of 25 mM L-histidine and 
various analogues (k obs for L-histidine = 0.11 min 1 ). 



relative fold 
cofactor & 0 bs discriminatioh 



L-histidine 


1 




L-histidinc methyl ester 


0.93 


1.1 


L-histidine benzyl ester 


0.76 


1.3 


ot -methyl-DL-histidine 


0.041 


24 


histidinamide 


0.025 


40 


glycyl-histidine 


0.006 


170 


histidino) 


0.003 


330 


3-methyl-L-histidine 


0.002 


500 


D-histidine 


0.001 


1000 


1-methyl-L-histidine 


< lO- 3 


> 1000 



The rate constant for HD2-promoted catalysis (k obs of 0.2 min* 1 , 50 mM 
histidine) is similar to that of natural self-cleaving ribozymes and corresponds to a rate 
5 enhancement of -10 million fold over the uncatalyzed reaction (k obs < 10' 8 min 1 under 
in vitro selection conditions). The dependence of the rate constant on histidine 
concentration is characteristic of the presence of a saturable binding site for histidine, 
although neither HD2 nor HD1 reach saturation even at 100 mM concentration of 
cofactor. The established specificity for particular cofactors, however, indicates that 
10 both catalysts do indeed form a histidine binding site. HD2 demonstrates greater 
activity with lower histidine concentrations, perhaps reflecting a greater binding 
affinity for histidine as would be expected due to its isolation from a low-histidine 
selection regiment. 

The pH-dependent activity profile for HD2 also implicates histidine as an 
15 integral component of the catalytic process (Figure 21b). The rate constant of HD2 is 
entirely independent of pH between the values 7 and 9. However, the activity of this 
enzyme drops precipitously at pH values that lie outside this optimum range. Most 
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revealing is the response of HD2 to low pH conditions. The k ohs values increase 
linearly with increasing pH between pH 4.5 and 5.5, giving a slope of approximately 
1 . This result is expected if the protonation state of a single functional group deter- 
mined the catalytic rate. Moreover, a rate constant that is half the maximum value is 

5 obtained at pH 6, where this chemical group will be half deprotonated. This value 
corresponds precisely with the pK, for the imidazole group of free histidine. Taken 
together, these results are consistent with a mechanism whereby the imidazole group 
serves as a general base catalyst for the deprotonation of the 2 '-hydroxy! group, 
thereby activating the oxygen for nucleophilic attack on the neighboring phosphorus 

10 atom. 

The loss of catalytic activity at higher pH values is not expected to be due 
to the protonation state of histidine, unless the pK^ of the imidazole group of a 
putative second histidine cofactor is dramatically shifted from its normal value. The (5- 
amino group of histidine, which has a pK, of greater than 9, conceivably could be 
15 involved in catalysis as well. However, it is expected to find a loss of activity with 
pH values in excess of 9 or less than 4.5 due to the significant level of deprotonation 
of T and G residues or protonation of C and A residues, respectively. 

Histidine was chosen as a candidate cofactor because of the potential for the 
imidazole side chain to function in both general acid and general base catalysis near 

20 neutral pH. This property is neither inherent to the four standard nucleotides of RNA 
nor to the remaining natural amino acids. As a consequence, histidine is one of the 
most- frequently used amino acids in the active sites of protein enzymes. For example, 
two active-site histidines are essential for the function of ribonuclease A from bovine 
pancreas, where both of these capacities are used to accelerate RNA cleavage. 

25 Although RNase A has long served as a model for the study of enzyme action, the 

specific roles that each active-site reside play in the catalytic process are still vigorous- 
ly debated (31). The classical view holds that the histidine at position 12 acts as a 
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general base for the deprotonation of the 2' hydroxyl, while the histidine at position 
119 acts as a general acid and protonates the 5' oxyanion leaving group. Breslow and 
others (25,47) have proposed that the role for histidine 119 instead may be to pro- 
tonate the phosphorane intermediate, thereby implicating general acid catalysis by the 
5 imidazole group as a priority step during the catalytic process. The data described 
herein indicate that the histidine cofactor for class II deoxyribozymes is not involved 
in a protonation step, but is functioning exclusively as a general base catalyst. 

In comparison to proteins, the more repetitive nature of monomelic units 
that make up nucleic acids limits both the formation of fine structure in folded 
10 polynucleotides and the chemical reactivity of RNA and DNA. The fact that a nucleic 
acid enzyme can co-opt one of the favorite chemical units of protein-based enzymes 
supports the notion that RNA could rally its limited structure-forming potential and, 
using the catalytic tools of modern protein enzymes, could produce and maintain a 
complex metabolic state. 

15 Materials and Methods 

In vitro selection and reselection 

In vitro selection was carried out essentially as described previously (5,7,47). 
The initial DNA pool was prepared by PCR amplification of the template 5-CTAAT- 
ACGACTCACTATAGGAAGAGATGGCGACATCTC(N) 40 GTGAGGTTGGTGTGGTTG 

20 (SEQ ID NOs: 53 and 54) (50 pmoles; N an equal probability of occurrence of the four nucleotides) in 
a 500-nL PCR reaction containing 400 pmoles of primer B2, 5 -biotin-GAATTCTAATACGACTCACT- 
ATrA (SEQ ID NO: 55), and 400 pmoles of primer 1, 5'-CAACCACACCAACCTCAC (SEQ ID NO: 
56), with 4 thermocycles of 94 9 C (15 sec), 50°C (30 sec), and 72X (30 sec). PCR reaction mixture was 
prepared as described previously (16). Amplified DNA was precipitated with ethanol, resuspended in 

25 binding buffer (50 mM HEPES (pH 7.5 at 23°C), 0.5 M NaCl, 0.5 M KC1, and 0.5 mM EDTA), and 

the solution was passed through a streptavidin-derivatized affinity matrix to generate immobilized single- 
stranded DNA 15 . The matrix displaying the pool DNA was repeatedly washed with 
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binding buffer (1.5 mL over 30 min), and subsequently eluted over the course of 1 hr 
with three 20-^iL aliquots of reaction buffer in which HEPES was replaced with 50 
mM histidine (pH 7.5, 23°C). In rounds 8-11, reaction time was reduced to 25-15 min 
to favor those molecules that cleave more efficiently. Selected DNAs were preciptita- 
5 ted with ethanol and amplified by PCR using primer 1 and primer 2, 5'-GAATTCTA- 
ATACGACTCACTATAGGAAGAGATGGCGAC (SEQ ID NO: 57), and the resulting 
PCR was reamplified as described above to reintroduce the biotin and embedded 
ribonucleotide moieties. 

Reselection of the class II deoxyribozyme was initiated with a pool of 10 13 

10 DNAs, each carrying a 39-nucleotide core that had been mutagenized with a degenera- 
cy of 0.21 per position. Similarly, HD2 reselection was conducted with an initial pool 
in which 26 nucleotides was mutagenized to a degeneracy of 0.33 per position. 
Individual from the final selected pools were analyzed by cloning and sequencing. 
The DNA pools were prepared for this process by PCR amplification using primer 2 in 

15 place of primer B2. DNA populations and individual precursor DNAs were prepared 
for assays as described previously (7). 

Deoxyribozyme Catalysis Assays 

All catalytic assays were conducted in the presence of 0.5 M NaCl, 0.5 M 
KC1, 0.5 mM EDTA. Single turn-over assays contained a trace amount (-50 nM) 

20 substrate oligonucleotide and an excess (1-10 jaM) DNA catalyst as described for each 
assay. The cofactor used was L-histidine unless otherwise stated. Reactions were 
terminated by addition to an equal volume of a solution containing 95% formamide, 
0.05% xylene cyanol, and 0.05% bromophenyl blue and stored on ice prior to gel 
electrophoresis. Termination buffers containing both urea and EDTA were incapable 

25 of completely terminating deoxyribozyme activity. 

Caged histidine experiments were conducted with intact dipeptides or with a 
concentration of hydrolyzed dipeptide products. Hydrolysis of dipeptides was achieved 
by incubating solutions containing 100 mM dipeptide and 6 N HC1 in a sealed tube at 
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1 15°C for 23 hr. Samples were evaporated in vacuo, coevaporated with deionized 
water, and the resuspended samples were adjusted to neutral pH prior to use. 

Catalytic rate constants (£ obs ) either were determined by determining the 
initial velocity of the reaction (16) or by plotting the natural log of the fraction 

5 substrate remaining over time, where the negative slope of the line obtained over 

several half lives represents fc obs . The uncatalyzed fate was determined by incubating a 
trace amount of 5' 32 P-labeled substrate under reaction conditions in the absence of 
deoxyribozyme at 23°C or at -20°C for 21 days. Comparative analysis of RNA 
phosphoester cleavage indicates that the rate constant for uncatalyzed RNA cleavage in 

10 the presence of histidine does not exceed the speed of substrate degradation due to 
radiolysis. It is expected that the maximum uncatalyzed rate for cleavage of the 
embedded RNA linkage does not exceed 10' 8 min'. This value is -10-fold lower than 
the value obtained in the presence of 1 mM Mg 2+ (7). 

The above description is for the purpose of teaching the person of ordinary 
15 skill in the art how to practice the present invention, and it is not intended to detail all 
those obvious modifications and variations of it which will become apparent to the 
skilled worker upon reading the description. It is intended, however, that all such 
obvious modifications and variations be included within the scope of the present 
invention, which is defined by the following claims. The claims are intended to cover 
20 the claimed components and steps in any sequence which is effective to meet the 
objectives there intended, unless the context specifically indicates the contrary. 
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